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ABSTRACT 

An experimental study was conducted to investigate the use of 
advanced composite materials in unstlffened thin-walled cylindrical 
structures subjected to hydrostatic pressure. Specifically, graphite/epoxy 
was looked at for potential application for small submersible pressure 
hulls. The investigation concentrated on the fabrication, instrumentation, 
testing, and evaluation of several 12-ply graphlte/epoxy cylinders of four 
different laminate configurations. 

The test specimens were 18 inches long with radii of 2.76 inches, and 
each was instrumented with 8 circumferentially and 2 axially mounted strain 
gages. Data collected from the strain gages throughout the tests were used 
in a modified Southwell type analysis and an analysis to deduce the 
buckling mode shape. The Southwell analysis method was shown to have a 
valid applicability to the buckling of laminated fiber-reinforced composite 
pressure vessels. A method was developed to use the bending strain data at 
locations around the centerline of the specimen to deduce the number of 
circumferential buckling waves. 

In addition, 1/2 inch ring specimen were cut from the excess ends of 
each specimen and tested under two-point compression, yielding information 
about the material properties of the laminates as fabricated. 

The hydrostatic tests were designed to verify the results of an 
analytical and numerical study that was conducted to determine the optimal 
laminate stacking sequences for buckling considerations. The test results 
indicate that the classical solution method presented by R.M. Jones 
produced accurate predictions of buckling strengths and correct relative 
rankings of the different laminate types. 
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NOMENCLATURE 



The following is a list of the special symbols and nomenclature used in 
this report. 



N. , N, 



y’ 



Njy = force resultants per unit length [Ibs/in] 

Ajj = laminate extensional stiffness matrix [Ibs/in] 
Bjj = laminate extension/bending coupling matrix [lbs] 
D^j = laminate bending stiffness matrix [lbs in] 
e = strain [in/in] 

K = curvature [1/in] 

= moment resultants per unit length [lbs] 
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X, 

U, 



0 = stress [Ibs/in ] 

E = equivalent modulus [Ibs/in^] 

R = mean radius [in] 
t = wall thickness [in] 

E = Young's Modulus [Ibs/in^] 

V = Poisson's ratio 

p = hydrostatic pressure (positive external) [psi] 
w = radial deflection (positive inwards) [in] 

L = length of the cylinder [in] 
n = number of circumferential buckling waves 
m = number of axial buckling waves 
y, z = axial, circumferential and radial coordinates 
V, w = axial, circumferential and radial displacements 
CLT = Classical Laminate Analysis 
GRP = glass reinforced plastic 
GFRP = graphite fiber reinforced plastic (also CFRP) 
W/D = weight to displacement ratio 
AUV = autonomous underwater vehicle (also UUV) 

GTE = coefficient of thermal expansion [in/in^F] 

FEA = finite element analysis 
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1. INTRODUCTION 



1.1. Potential Benefits of Composite Pressure Hulls 

Since the 1960 's fiber reinforced composite materials 
have been proposed for the construction of submersible 
pressure hulls. The potential benefits that can be offered by 
composite materials over metals include: 
o structural weight reductions; 
o fabrication of novel hull geometries; 
o incorporation of acoustic silencing; 
o electromagnetic signature reduction; 
o virtual elimination of corrosion. 

The purpose of this project is to gain insight into the use of 
graphite/epoxy for this application, through the systematic 
testing to hydrostatic pressure collapse and subsequent 
analysis of several thin-walled laminated composite cylinders, 
of various laminate types. 

Weight Reduction. Composite materials offer the 
opportunity to make reductions in the portion of a 
submersible's weight that has to be devoted to the hull 
structure. In fact, submersible designers use as a measure of 
hull efficiency the ratio of the hull weight to the weight of 
the water the hull displaces. This ratio is referred to as the 
"W/D" . It must be less than one for the submersible to float, 
and a smaller "W/D" means that the hull is more efficient. 



9 




Figure 1.1. Relationship of weight to depth ratio to collapse depth 
for several material systems. 

because for a given depth it can either carry more payload or 
use less energy. Values approaching 0.30 are required for the 
design of small vehicles with adequate allowances for mission 
capabilities. Figure 1.1 shows graphically the relationship of 
W/D to collapse depth for several materials [1]. From this 
figure it is seen that composite materials can offer 
substantial increases in collapse depth at any given W/D, or 
that for a specified collapse depth requirement composites 
enable large reductions in structural weight [2]. 

The major structural weight savings obtained via use of 
composite material construction can be directly utilized in 
other mission enhancing areas. Table 1.1 compares the typical 
weight allocations for an Autonomous Underwater Vehicle (AUV) 
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[1] with current U.S. Navy attack submarine (SSN) and fleet 
ballistic missile submarine (SSBN) design practices. Small 
submersibles are very volume critical, and designers always 
have more and larger payloads than they are capable to carry. 
Weight reductions in the hull structure can be translated into 
increased payload weight and volume capacity. For the many 
small submersibles where speed or range increases would be 
more beneficial, portions of the structural weight savings 
could be channeled into larger propulsion machinery or fuel 
carrying capacity. For submersibles that require the 
capability to reach greater safe operating depths, the 
incorporation of composite materials could make the hull 
structure stronger without giving up established payload, 
propulsion and fuel allotments. 



CATEGORY 


AUV 


SSN 


SSBN 


1. 


Hull Structure 


30.6 


34 


32 


2. 


Variable Ballast 


3.7 


12 


12 


3. 


Propulsion 


32.3 


25 


13 


4. 


Payload 


13.4 


2 


17 


5. 


Fixed Ballast 


1.5 


5 


12 


6. 


All Others 


18.5 


22 


14 



Table 1.1. Typical weight allocations of submersibles. 

Novel Hull Shapes. The methods of fabricating composite 
material hull structures make it easier for the designer to 
create complex geometrical shapes and incorporate advanced 
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sandwich and double hull structural configurations. The use of 
novel structural geometries such as figure-eight and swedge 
hull shapes, allow for improved structural efficiency by 
optimizing the usage of space and volume as compared to 
traditional cylindrical and spherical hullforms [3]. 

Filament winding fabrication techniques allow for complex 
streamlined axisymmetric shapes to be constructed with 
relative ease. 

Signature Control . A major advantage of laminated 
composite pressure hulls is in the reduction of acoustic and 
non-acoustic signatures. This can be crucial in missions such 
as mine hunting and covert reconnaissance. Composites have 
seen years of usage as sonar windows due to their inherent 
acoustic transparency [4]. Also, by innovative utilization 
of low density sandwich cores and other acoustic damping and 
decoupling lamina and anechoic coatings, internal noises from 
machinery and electronics can be silenced to provide improved 
stealth. Because they are non-metal lie, these materials are 
also non-magnetic and electrically non-conductive, greatly 
lowering a submersible's susceptibility to enemy weapons and 
sensors. 

Corrosion Control . Corrosion, which is a major problem 
for metals in the marine environment, can be virtually 
eliminated with properly selected composite material systems. 
However, in completeness, it must be mentioned that initial 
attempts at utilizing composites, specifically glass 
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reinforced plastic (GRP), found that seals around hull 
penetrations, hatches and viewports were susceptible under 
cyclic loading. The high deformation of the GRP as compared to 
metallic inserts left cut fibers vulnerable to exposure to 
water at high pressure. Because of this designs have been 
developed where cylindrical composite pressure hulls are used 
in conjunction with specially designed titanium end closures. 
In these designs all penetrations are made through the 
titanium end closures thus avoiding the difficult problem of 
penetrating the laminated composite pressure hull. Due to the 
weight of the metal end closures and because the unstiffened 
hulls are buckling critical, the overall hull weight may not 
show as significant reductions as predicted for all-composite 
hulls. Notwithstanding this, all the benefits presented by 
composite materials and the fact that any weight reduction is 
precious to a submersible designer, these advanced materials 
are well suited for utilization in pressure vessel design and 
construction. Graphite fiber reinforced plastics (GFRP) 
because of their higher modulus, lighter weight, and 
extraordinary fatigue properties generally offer more long 
term benefits than GRP. 
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1.2. Purpose of Hydrostatic Testing 



Four areas of major concern in the design of laminated 
composite pressure hulls are: 

o fabrication procedures, 
o repairability and serviceability, 
o inspection and certification, 
o in-service hull monitoring. 

The present investigation, consisting mainly of the 
hydrostatic testing of 12 -ply graphite/epoxy cylindrical 
specimens of various layups, is aimed at providing data and 
techniques to better facilitate the addressing of these 
concerns . 

Prior to and concurrent with this test program, an 
analysis, both classical and numerical, was performed to 
predict the buckling characteristics of the graphite/epoxy 
cylinders that were tested here [ 5 ]. In this analysis, 
closed form solutions were able to be used for the cylinders 
constructed of orthotropic laminates. These laminates do not 
exhibit any coupling between extension and shear or bending 
and twisting stiffnesses (i.e. they have no Ai 6 , A26/ Bu, B26, 
Di 6 , or D26 terms in their ABD matrices, which relate stress to 
strain) . Finite element analysis was also conducted on these 
orthotropic laminates and additionally on laminates with angle 
plies. These angle-ply laminates generally have some coupling 
between extension and shear or bending and twisting due to the 
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presence of these additional stiffness terms. 

Each specimen was instrumented with eight equally spaced 
circumferential strain gages. The average of these strains was 
compared to the axisymmetric strain that is predicted by 
theory and was also used as the membrane strain for a 
Southwell -type analyses to predict perfect shell buckling 
characteristics from experimental data from imperfect 
specimens . 

In summary, the main work on this thesis was to conduct a 
large series of experiments in an effort to verify the results 
of theoretical elastic stability analyses methods that have 
been adapted for laminated cylindrical shells subjected to 
hydrostatic pressure. Since shell buckling characteristics are 
highly susceptible to leiminate imperfections, geometrical 
deviations, and residual stresses, the results of this testing 
depended greatly on the effectiveness of the specimen 
manufacturing method. The fabrication technique chosen for 
this project consistently produced specimens of high 
circularity, but the following defects were found to varying 
extents in all the cylinders: wrinkled fibers, non-uniform 
resin distribution, voids and del ami nations. In several of the 
specimens noticeable defects, especially severely wrinkled 
plies, were apparently responsible for large degradations in 
the buckling strengths of these cylinders. 

Although these circular cylindrical shells exhibit 
significant buckling sensitivity to geometric shape 
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imperfections just like isotropic shells, in general, 
laminated composites offer the opportunity to reduce these 
effects by selecting appropriate fiber orientations. An 
overall goal of this entire analytical and experimental test 
program is to develop and verify a method to optimize the 
selection of ply orientations in order to achieve the most 
efficient hull for specific requirements. Before optimization 
routines can be applied to numerical or closed form solutions, 
there must be confidence that these solutions are fairly 
accurate and can at least rank various laminate configurations 
correctly. Four different laminates were chosen to be 
evaluated by experimental test methods in an effort to test 
the analyses; 

o 1:1 specially orthotropic - [QOa/Oajs / 
o 2:1 specially orthotropic - [904/02]s , 
o quasi -isotropic - [902/(±30)2]s / 

o interleaved (unsymmetric ) orthotropic - [902/0]4 t . 
By testing multiple specimens of each type, and comparing the 
rankings with the results of the analysis, insight into the 
correctness of the theoretical methods was obtained. 

Another significant contribution of this hydrostatic 
testing results from the extensive strain gage instrumentation 
that was utilized on each cylindrical specimen. On the 
centerline of every specimen were placed eight circumferential 
strain gages at 45° separation and two axial strain gages at 
180° separation. Readings from these gages provided the 
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ability to validate the assumptions of buckling prediction 
theories, to investigate the applicability of alternative 
methods of stability determination, and to provide insight 
into the mode shapes into which the cylinders buckled. 

In the 1930 's, R.V. Southwell proposed a method to 
determine the buckling strength of an ideally straight column 
by analysis of test data taken from a column with initial 
curvature [6]. In the 1950 's, G.D. Galletly and T.E. 

Reynolds proposed that a extension of this method may be used 
to determine the elastic stability of ring-stiffened cylinders 
[7]. The applicability of this method to determine the 
stability of composite Icuninated cylinders was investigated 
using the strain data collected during the hydrostatic 
pressure tests. 

By isolating the bending strain from the membrane strain 
at the eight locations around the centerline of the specimens 
and applying elementary concepts that relate the local radial 
deflections to the changes in circiomference, a deduction was 
made of the buckling mode shape of each of the specimens. 

These mode shape deductions are compared to those assumed and 
predicted by the analysis methods. 
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2. ANALYTICAL BACKGROUND 



This section addresses the significant issues in material 
science and in structural stability pertinent to the failure 
of laminated, fiber reinforced composite cylinders subjected 
to external hydrostatic pressure. Designing with laminated 
composite materials offers the opportunity to achieve specific 
structural properties through proper selection of: 

o the system of resin, fibers and enhancing agents; 
o processing and production techniques; 
o individual ply orientations. 

To make these selections the designer must consider concerns 
ranging from the micromechanical bonding of the fibers and the 
matrix, to the ability of the structural fabrication technique 
produce highly circular cylinders. This chapter attempts to 
bring together some of the information and experience that 
have been gained in these broad fields. Only by making the 
material scientist and the naval architect aware of the 
other's knowledge and concerns, can all the benefits presented 
by the utilization of composite materials in submersible 
pressure hull design and construction be realized. 

The discussion first addresses the material science 
concerns, beginning with the micromechanical factors 
controlling the compressive failure of fiber reinforced 
laminated composites, and leading to a macromechanical 
treatment of Classical Laminate Theory. 
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The discussion then covers the major structural concerns, 
starting with the axisyimnetric compression of circular 
cylinders, progressing through the compression of thin rings 
to the elastic stability theories for long cylindrical shells. 

The two fields merge with a coverage of a classical 
buckling analysis method for laminated cylinders, that is 
based on Donnell's stability theory. Figure 2.1 shows the 
coordinate system and dimensional nomenclature used in this 
analysis. Next, brief coverage is given of some recent finite 
element studies on this subject. 

Finally, a modified Southwell -type method is presented 
that should be universally applicable to a wide range of 
isotropic and anisotropic material systems. Additionally, a 
method is proposed whereby the bending strain data can be used 
to deduce the nximber of circumferential waves in the actual 
buckling mode shape of a cylindrical specimen. 




Figure 2.1. Diagram of the coordinate system used In the analysis. 
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2.1. Compressive Strength of Composite Laminates 



The compressive response of fiber-reinforced composite 
materials is still one of the least understood in the field of 
composites today. The compression failures of composite 
laminates can be categorized into the following mechanisms: 
"global Euler buckling, microbuckling, transverse tension, 
fiber kinking, fiber compression failures, matrix compression 
failures, and delamination." [8] The compression failure 
modes of fibers, matrices and unidirectional composites are 
illustrated in Figures 2.2 and 2.3, which were taken from 
reference 11. The predominant mechanism is determined by 
factors such as the individual properties of each of the 
constituent materials, the orientation of the plies within the 
laminate, the geometry of the specimen, the method of load 
application, any waviness and misalignment of the fibers, the 
presence of voids and discontinuities in the matrix, and any 
stress concentrations. 

If the fibers are weak enough in compression, they could 
fail due to compression before the onset of buckling and cause 
ultimate failure of the laminate. But unlike tension, the 
compressive strength of most composites is limited more by the 
properties of the matrix than the fiber reinforcement. 
Microbuckling can begin to occur because the resin is soft, 
initiating ultimate failure, even when the fibers are not weak 
in compression. Other reasons the matrix may not be able to 
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provide adequate lateral support to the fiber are: proximity 
to a free surface, the formation of lateral splitting, a high 
void content, or environmental softening. Desirable matrix 
properties include a moderate amount of stiffness and 
toughness, but not an excess. Elastic microbuckling is the 
dominant failure mode when either the fibers have no initial 
misalignment or the matrix is very soft. The low compressive 
strengths of fiber reinforced composites may also be the 
result of kink band formation, an occurrence that was 
predicted by theory and is now possible to observe 
experimentally. This type of compression failure results from 
fiber misalignment and involves plastic deformation of the 
resin. 

Microbuckling. Fiber buckling theories are based on the 
stability of fibers in a flexible matrix. The fibers may 
buckle in either of two modes: (1) when the fiber volume 

fraction is greater than 30 percent, the fibers tend to buckle 
in phase (shear mode), and (2) when the fiber volume is less 
than 30 percent, the fibers tend to buckle out-of-phase 
(extensional mode) [9]. In a study of the effects of 
various constituent properties, fiber arrays, and bowed, 
unbonded and misaligned fibers, Greszczuk concluded that the 
microbuckling models cannot predict the failure of non- 
ideal ized, actual composites. He stated further that low- 
modulus resins result in microbuckling failure, but when 
higher modulus resins are used the failure mode changes and 
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the compressive strength of the fiber reinforcement controls 
the failure [10]. 

Kink Band Formation. When the resin system is adequately 
stiff and there exists any fiber waviness, regardless of the 
mode of incipient failure, final failure occurs as the 
formation of kink bands. Chaudhuri conducted an analytical 
study to identify the basic failure mechanisms of 
graphite/epoxy cylinders with thick walls under hydrostatic 
compression [11]. He identified the dominant compressive 
failure mode for thick composite cylinders to be the 
"formation and propagation of fiber kink bands at the 
microscopic level, triggered by the fiber misalignment defects 
formed during the manufacturing process, leading to the shear 
crippling failure at the macroscopic level." He found that 
this kinking was the lowest energy compression failure mode at 
the microlevel of graphite/epoxy composites, and was limited 
by the initial fiber misalignment, the ultimate fiber strain 
and the two transverse shear moduli, Glt and Git, of the 
laminate . 

Chaudhuri ' s study is important because it shows that kink 
band formation is a distinct and lower energy failure mode 
than elastic microbuck ling, not due solely to the effect of 
plasticity in the presence of initial fiber misalignment. 
Chaudhuri also noted that the ratio of hoop to axial layers 
did not strongly influence the compressive strength of the 
thick-walled (R/t<6) composite cylinders he investigated. 
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Higher toughness usually means fiber kinking will occur 
at lower strains, but lower toughness (i.e. higher modulus) 
usually means that there is less resistance to delamination 
propagation. Therefore, there must be a tradeoff between the 
delamination associated with higher modulus resins and the 
earlier (i.e. at lower strains) microbuckl ing initiation found 
with tougher resins. 

Shuart presented a linear analysis of short-wavelength 
modes shapes for off-axis and angle-ply graphite/epoxy 
laminates and a non-linear analysis of laminate stresses, out- 
of-plane displacements, and interlaminar shear strains 
[12]. He showed that significant interlaminar shearing 
strains develop before short-wavelength buckling, and he 
described a compression failure criterion that included both 
these failure mechanisms. 

Even though fibers do not creep, all polymer matrices 
exhibit viscoelastic creep behavior which is affected by water 
absorption. Mechanical impact damage and exposure to high 
temperatures can aggravate this problem. Irion and Adams 
presented results that indicate that longitudinally reinforced 
graphite composites, in general, do not creep any significant 
amount. They found that transversely reinforced graphite 
composites exhibited long periods at high strains with no 
creep followed by abrupt jumps. Because few of their specimens 
attained steady-state creep and also due to the fact that the 
strain versus time plots displayed extensive amounts of 
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irregularities, they inferred that local failures were 
occurring almost continually, but after initiation the 
failures were quickly arrested [13]. 

Transverse tensile stress in unidirectional composites 
has been proposed as a reason experimental compression 
failures occur at much lower loads then predicted by theory. 
Even though this stress is small, unidirectional composites 
have low transverse strength. 

In view of all this the expected failure sequence for 
laminated shells constructed from plies of unidirectional 
graphite fiber-reinforced composites with fiber volume 
fractions greater than 30%, and loaded in hydrostatic 
compression can be summarized as follows: "microbuckling in a 
shear mode (in-phase buckling) will occur, initiated by 
regions of local inhomogeneities (voids, stress 
concentrations, areas of weak matrix, free boundaries), and 
lead to final failure in the form of fiber kinking or fiber 
overstressing in bending." [8] 

A major limiting factor on the compression strength of 
laminated composites is the technique of fabrication. 
Manufacture of the specimen has significant effect on the 
composite's compressive response because the presence of local 
inhomogeneities and defects, which are difficult to model and 
otherwise account for, influence the compressive failure 
mechanisms on the microstructural level more than any of the 
other stress state failure mechanisms. Therefore, the chosen 
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method of fabrication must result in a low void content (less 
than 2 %), because otherwise the fibers located adjacent to a 
void will not have sufficient lateral support. Laminates must 
also be fabricated with very low amounts of fiber misalignment 
(compressive strength can be highly sensitive to fibers out of 
alignment as little as 2°). Some research suggests that 
processing with high fiber content (greater than 60%) may 
result in lower void contents and, therefore, more precise 
fiber alignment. 

The 5208 epoxy resin used in the fabrication of the 
cylinders for this present study has a high modulus, with 
failure strains of 1.4%. This low resin toughness has a 
detrimental effect on the progression of the failure. Studies 
have shown that compression failure of brittle laminates, such 
as T300/5208, is quite sudden, and that it is very difficult 
to arrest the resulting fiber kinking [14]. Sinclair and 
Chamis conducted compression tests on unidirectional T300/5208 
graphite/epoxy laminates and their results were divided into a 
high-strength group where failure was predicted by fiber 
compressive failure, a medium strength group that was 
predicted by flexure or del ami nation, and a low strength group 
that was predicted by Euler Buckling [15]. 
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2.2. Classical Laminate Theory 



Before discussing the elastic stability of composite 
pressure hulls, a brief overview of Classical Laminate Theory 
(CLT) is given. The concept of equivalent moduli for 
orthotropic cylinders is also developed to enable 
relationships designed for isotropic rings and cylindrical 
shells to be adapted, without excessive complexity, to the 
case of the laminated composite rings and cylindrical 
specimens tested in this research. 

For laminated orthotropic materials the constitutive 
relationship between stress and strain is complicated by the 
directional nature of the individual plies in the laminate. 
The elastic modulus for isotropic materials is replaced by an 
ABD matrix that contains information on the extension, shear 
and bending behavior of the laminate and the coupling between 
the extension, shear and bending. The normal presentation of 
this constitutive relation is 
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where 



w 



^ij ^ ^k-1 ^ 



( 2 ) 




(3) 




(4) 



In general, the Bij terms couple the bending and stretching 
effects. If the laminate is symmetric about the middle plane, 
there is no such coupling and thus Bij=0. In specially 
orthotropic laminates, which are cross-ply Iciminates composed 
of plies oriented at 90° relative to each other, there is no 
coupling between bending and twisting effects, and, therefore, 
Ai6=A26— Di6 = D26 = 0 . 

To enable comparison of the relationships in the 
following sections for the axisymmetric compression and 
elastic stability of long composite cylinders with the 
experimental results, equivalent properties for an orthotropic 
cylindrical shell were developed. 

The membrane constitutive relations for an orthotropic 
material, where the bending response and extension/shear 
coupling are neglected, are 
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This is inverted to solve for the strains, such that 
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For the case of hydrostatic loading. 
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and thus 
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and the equivalent modulus in the hoop direction is therefore. 
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Similarly, the equivalent modulus in the axial direction is 



■^ H -^22 -^12 
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(11) 



and the equivalent Poisson's ratio is 



V = 




= V 



12 



( 12 ) 



These expressions can be used later with the experimental 
results, to check the effective properties of the cylindrical 
specimens as fabricated. Note that these relations only apply 
to the case of hydrostatic loading, for other loading 
conditions other equivalent moduli must be developed. 



30 



2.3. Azi symmetric Compression of Cylindrical Shells 

It has been said that axisymmetric yield failure is the 
most important mode of failure of externally loaded pressure 
vessels because a "sensible design should eliminate failure 
due to instability." [16] This is lofty goal for thin- 
skinned pressure vessels because initial geometric 
imperfections on the order of a tenth of the thickness can 
cause a catastrophic loss of buckling resistance of pressure 
vessels when subjected to external water pressure. 

One of the earliest solutions for the axisymmetric 
deformation of circular cylinders under external loading was 
published by von Sanden and Gunther in 1920 [17]. Their 
original work unfortunately had two signs interchanged in the 
formula for circumferential stress, which can give incorrect 
values if used unknowingly. Also, these authors did not take 
into account the beam-column effect and its associated non- 
linear deformations, that was later introduced by Salerno and 
Pulos [18] and refined by Ross [23] into the form 




dK^ dx^ 
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where 



= 3(l-v^) 

02 = pR^ + 

2Et 12(1-v2) 

pi?2(l-^) 



( 14 ) 



Ross gives the complete solution, but here is presented only 
the solution for our case of a circular shell element, 
stiffened by equal -strength frames located at its ends 
(x=0,L), and deformed symmetrically about the midspan. 
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where 



i=\ = aV (1-a^p^) Fj = a\/ (1+a^p^) 

and A, and Aa are functions of Fi and Fa that reflect the 
characteristics of the end frame and the boundary conditions 
assumption that dw/dx=0 where the edge connects to the end 
frame (x=0,L). The variation of the hoop and axial stress 
along the shell is due to the radial contraction, the 
Poisson's effect of the axial compression, and the beam-column 
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effect. These are expressed in the following relationships for 
hoop stress and axial stress, respectively, where the last 
term is the beam-column effect and is positive in the inner 
fibers and negative in the outer fibers. 



HoopStress: Oy 
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AxialStress: 
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2.4. Elastic Stability of Long Cylindrical Shells 

The discussion of the elastic stability of laminated 
composite shells begins with a review of classical buckling 
theories for isotropic shells of revolution. This paper 
addresses only cylindrical shells whose thickness t is small 
compared to the other dimensions of the shell, especially the 
radius of curvature. The discussion begins with a thin ring 
under symmetric two-point compression and works toward the 
theories for a cylinder of length L and radius R. 

The theory for thin isotropic circular rings under 
compressive loading was presented by S.P. Timoshenko in his 
classic text on elastic stability [19]. A summary of that 
analysis is given, and extended for the case of laiminated 
composite rings. The relationships obtained here will be used 
in the analysis of the results obtained from the two-point 
compression testing of 1/2 inch ring specimens that were cut 
from the excess ends of each cylinder. 

Consider a isotropic circular ring of unit width, 
subjected to a two-point compression load P as shown in figure 
2.3. The bending moment at A (and at B) is denoted Mo and the 
moment at an arbitrary location 0 is 

M = + ■— (1-COS0) (19) 

The differential equation for the radial deflection of this 
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Figure 2.4. A ring of radius R compressed by two forces P acting 
along a diameter. 



ring at any point 0 due to this loading is 
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From the symmetry of the problem where 

-^ = 0 at 0 = 0 and 0 = -^ 
dd 2 



( 21 ) 



the general solution of equation 20 is found to be 
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Timoshenko then applied Cast igliano ' s theorem to the strain 
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energy expression to get 



Mo = 



^(2_i) 

2 n 



(23) 



Of interest to us is the radial deflection at the point 
of load application (0=n/2), which from equation 22 is 



w 
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(24) 



where w < 0 indicates the shortening of the diameter or 
deflection inward toward the center. 

For the case of a thin (R>>t) laminated composite ring, a 
two-point, symmetric compressive force has a negligible direct 
effect on the stress on the inner and outer surfaces at point 
A, and this stress, therefore, is almost entirely due to the 
bending moment resulting from that applied point force. 
Additionally the compressing ring is assumed to undergo 
cylindrical bending ( k = 0 ) and not anticlastic bending. 
Under these conditions, equation 1 can be simplified to 




With cylindrical bending and negligible Poisson's effect, we 
get that 
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From the definition of curvature as the second derivative of 
the displacement we see that effectively. 
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(27) 



Substituting this effective bending stiffness of the laminated 
ring, Dja’b = El , into Timoshenko's formula for an isotropic 
ring, equation 24, and solving for D 22 gives 



= 0.0744 






(28) 



On the inner and outer surfaces the strain is related to the 
curvature, such that 






€ + K .Z 



(29) 



where is negligible for this loading, as explained above, 
and Zupper = t/2 on the surfaces. Combining equations 23, 26 
and 29 gives a relationship for the hoop strain on these 
surfaces in terms of the applied compressive load P, such that 



€ = tO.1817 — 

2 D,,b 



(30) 
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Solving for D 22 gives 



D 22 = TO. 0908 



b 




(31) 



Therefore, we now have a way to approximate the D 22 of the ring 
with the slope of the load versus displacement plot, equation 
28, and with the slope of the load versus strain plot, 
equation 31. These experimentally determined D 22 ' s can be 
compared to the one obtained from CLT and the properties of 
T300/5208. The complete ABD matrices for each of the laminate 
types were developed in this way and are contained in Appendix 
E. This comparison provides insight into the actual bending 
stiffness properties of each of the laminate types as 
fabricated . 

Next, consider now how the radial deflection in the ring 
under the two-point, symmetric loading would be changed if a 
hydrostatic pressure of magnitude p was applied. This uniform 
external pressure produces a compressive force S = pR in the 
ring of unit width. This force combined with the symmetric, 
two-point loading considered above produces the deflection 
expressed by 
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This means that each term in w increases as the compressive 
force S increases, becoming indefinite when 



S 



(n^-1) El 



(33) 



In this case the critical compressive force occurs when n=2, 
therefore 



2pR^ cos26 

9nEI (34) 

^ciit 



This thin ring analysis can now be extended to the 
buckling of long circular cylinders under uniform external 
pressure. The critical pressure Pcrit on an elemental ring of 
unit width of the cylinder is obtained by considering 
poisson's effects that were not a part of the isolated ring 
above, i.e. substitute E/(l-vM for E and t^/12 for I. 



Pczit 
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(35) 



This calculated critical pressure is applicable as long as the 
compressive stress does not exceed the material's proportional 
1 imit . 



One of the most significant contributions to this subject 
was the general theory of shell instability developed by W. 
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Koiter in 1945 [20]. In his thesis, Koiter demonstrated 
conclusively that initial geometric imperfections play a 
dominant role in reducing the buckling load of shell 
structures. This helped explain the scatter in experimental 
buckling test results and the general inability of researchers 
to experimentally attain theoretically predicted critical 
pressures. The greatest detriment to achieving theoretical 
buckling pressure predictions is imperfections in the shape of 
the shell. Koiter showed that in isotropic shells an 
imperfection amplitude of just 10% of the wall thickness of 
the shell can result in as large as a 40% reduction in its 
buckling strength. This can be uncontrollable in thin-skinned 
shells with R/t>>10 like those studied here. More recent work 
by Koiter and others has led to theories where the effect of 
small but finite imperfections are included in the stability 
analysis. These give the designer the ability to account for 
geometric imperfections that can be assumed to be present in 
the as-manufactured pressure vessel. Noteworthy of these is an 
equation derived by Kendrick [21] which describes the 
stress in an isotropic cylinder due to eccentricities, such 
that 
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where Co is the amplitude of the initial deviation in a mode 
shape cos(n0). The most important term in this equation is 
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( 37 ) 



Co ( 



Peri t P 



which shows that the hoop stress is directly proportional to 
the amplitude of the initial eccentricities in the cylinder. 

It also verifies that a failure envelope exists for a cylinder 
under external hydrostatic pressure because the stress becomes 
singular as the critical buckling pressure is approached. 
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2.5. Classical Buckling Analysis of Composite Shells 

This section reviews some of the more prominent theories 
on the buckling stability of laminated composite cylindrical 
shells. These theories are referred to as "classical" because 
they all relate the concepts of Classical Laminate Theory 
(CLT) to one or more of the classic elastic stability 
theories. The analysis begins with circular cylinders 
constructed of orthotropic layers and is based on work 
presented by R.M. Jones. Next, the analysis considers the 
introduction of angle plies and is based on work by Booten and 
Tennyson. Lastly, the analysis proposes the application of 
optimization techniques to laminate orientation selection and 
is based on the recent work of several different authors, most 
notably Sun. Please note that subscripts following a comma 
represent partial differentiation by that subscript. For 
example, the term N,,, is equivalent to the expression 3N./3x. 

In 1968, Robert M. Jones presented an exact solution for 
the buckling of a circular cylindrical shell with multiple 
orthotropic layers under any combination of axial compression 
and lateral pressure [22]. Jones showed that previous 
attempts to predict the buckling behavior of multilayered 
shells were erratic because they neglected the coupling 
between bending and extension, which his theory considered. He 
based his work on Donnell's stability differential equations 
for circular cylindrical shells, with an implied membrane 
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prebuckled shape and simply supported end conditions. 

By applying the Ki rchhof f-Love hypothesis, the three- 
dimensional orthotropic elasticity relations were reduced to 
two-dimensional stress-strain relations for the layer. The 
variations in stresses in the k"'' layer during buckling are 



where 6i, 62, and €3 are the variations in reference surface 
strains, and ki, K2, and ks are the variations in reference 
surface curvature, such that 
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and the variations in strains during buckling are 
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and 
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Thus the variations in stresses in the k"" layer can be 



expressed 
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io* = On (Ci + ZKi) + On (€2 + 2X2) 



= O 12 (€^+ZK^) + O 22 (€2 + 2 X 2 ) 



( 42 ) 



= Osi (€j+ZKj) 

These variations of stresses over the shell layers during 
buckling are then integrated to give the variations of forces 
and moments during buckling; 
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These force and moment resultants, which account for the 
coupling between bending and extension, only require four 
independent material properties, E,,", Eyy“, G.y‘‘, and Viy", per 
layer. 

L.H. Donnell proposed the following stability 
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differential equations for circular shells subjected to 
combined axial compression and lateral pressure [23], 




( 49 ) 
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Substituting into these equations, the variations of the 
in-plane forces per unit length and the variations of the 
moments per unit length, equations 43-48, together with the 
previous stress-strain and strain-displacement relationships, 
equations 40-42, gives three stability differential equations 
in terms of the displacements, u, v and w: 



where Li, are appropriate partial differential operators, see 
reference 22. 

There are sixteen possible force and geometric boundary 
conditions of which a combination of four constitutes a 
defined problem. Our main interest is in simply supported 
ends, such that. 
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These four conditions are satisfied by the following buckling 
displacements: 



u = u cos 
V = V sin 
w = w sin 




(53) 



where the overline indicates the amplitude. Placing these 
expressions into equation 51, leaves the following 
relationship, whose determinant must be zeroto obtaino be 
nontrivial solution: 




(54) 



This results in the following stability criterion, which we 
shall refer to as "Jones' Equation": 
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where 
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The buckling load under hydrostatic pressure is obtained by 
substituting the following load factors: 



N = ( 62 ) 

y t ^ 2 t 
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Jones' contribution was instrumental because it 



introduced the effect of the coupling between bending and 
extension that exists due to the presence of the multiple 
orthotropic layers in laminated, composite shells. The more 
difficult questions of initial geometrical and material 
imperfections and prebuckling deformations were not 
considered . 

In 1978, M. Booten and R.C. Tennyson conducted a study 
that considered the behavior of laminated anisotropic 
cylinders with small axi symmetric shape imperfections 
[24], Their analysis used exact boundary conditions and 
nonlinear prebuckling deformations, and comparisons were made 
with experimental data. These authors concluded that the 
effect of varying the fiber angle on the interactive buckling 
behavior of laminated cylinders with "perfect" geometry and 
"less than perfect" axisymmetry is significant. In a later 
work with Hansen [25], Tennyson demonstrated that the 
buckling strength of an axially loaded cylindrical shell could 
be increased significantly through a judicious choice of 
laminate configuration. This has naturally lead to the topic 
of optimization. 

Composite materials offer the designer the ability to 
tailor material properties to the specific application. In 
filamentary composites, such as graphite-epoxy, this involves 
ply angle selection, where designers aspire for the capability 
to select optimal ply angle distributions for specific 
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geometries and load conditions. In the past ten years several 
studies have been conducted into laminate optimization for 
cylindrical shells. 

In 1981, Y.S. Nshanian of the then U.S.S.R. and M. Pappas 
examined the optimization of the ply angle distribution across 
the thickness of composite shells for buckling and vibration 
considerations [26]. They proposed methods to determine 
the optimal ply angle variation through the thickness of 
symmetric angle-ply shells of uniform thickness. Using a 
mathematical programming algorithm, the procedure attempted to 
maximize the buckling load of a thin, simply-supported, 
circular, cylindrical, angle-ply shell. Their results showed 
significant gains using optimal variable ply angle 
configurations, as compared to ones with an optimal constant 
ply angle. But reference shells with variable but non-optimal 
angle-ply configurations (0°/90° and 90°/0°) were often near 
optimal in performance and always superior to optimal 
conventional angle-ply shells. 

In 1987, Sun and Hansen attempted to maximize the 
buckling load of a laminated-composite, circular-cylindrical 
shell as a function of ply orientation [27]. They used a 
two-step optimization routine with a first step random search 
to determine an initial guess and a second step systematic 
search based on work by Powell [28] . The optimization 
routines were applied to the finite difference solutions to 
von Karman-Donnel 1 non-linear shell strain-displacement 
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relations. Experiments were conducted on selected, optimal, 
four-ply, graphite/epoxy cylinders and compared to certain 
reference laminates. The experimentally determined buckling 
load under external pressure for the optimized cylinder was 
3.19 times as high as the [0/90]s reference cylinder, but no 
[90/0]s reference cylinders, with hoop plies on the outer and 
inner surfaces, were tested under external pressure. 

In 1989, Sun presented another optimization study using 
linear theory and the same two-step optimization method 
[29]. The study assumed the cylinders were constructed of 
N orthotropic plies with the same thickness, with the i'^*' ply 
containing an equal number of fibers in the +6i and -0t 
directions uniformly distributed through the ply's thickness. 
This eliminates the stiffness coefficients of the laminate Aie, 
A 26 / Bi6, B 26 , Di 6, and D 26 . This assumption is valid for 
laminates with many +0 and -0 plies of thicknesses much 
smaller than the thickness of the laminate. It serves only as 
an approximate result for laminates with relatively few plies, 
such as the 4-ply ones Sun tested and the 12-ply laminates 
tested here. 
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2.6. Finite Elements Analysis of Composite Shells 

In 1984, Abu-Farsakh and Lusher investigated the buckling 
behavior of glass-reinforced cylindrical shells under combined 
axial compression and external pressure [30]. Initial 
imperfections, load-deflection response, and buckling modes 
and loads were determined experimentally for 20 specimens and 
compared to the results of a finite element analysis (FEA). 

The FEA predicted lower buckling loads than were 
experimentally measured. They concluded that the size of the 
mesh must produce a higher nvimber of divisions 
circumferentially than the number of circumferential buckling 
waves to obtain reasonable results. 

As a prelude to the hydrostatic testing outlined in this 
report, extensive finite element analysis work was completed 
by T. Perry at the University of Massachusetts, Lowell, 
Massachusetts [5] . The suitability of three finite element 
software packages were evaluated for use in this project, EMRC 
NISA, Algor’s SUPERSAP, and HKS ' S ABAQUS. ABAQUS, which has 
anisotropic material capability for buckling and non-linear 
analysis, proved to be the best of the three packages. ABAQUS 
finite element models for specially orthotropic, orthotropic, 
quasi-isotropic, and isotropic cylinders were run using the 
material properties of T300/5208 graphite/epoxy and simply 
supported boundary conditions. The models were created using 
16 circumferential by 10 longitudinal (16cxl01) parabolic 5 
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degree of freedom/node elements for a total of 2560 degrees of 
freedom (dof). The FEA gave stiffer results by as much as 25% 
over the classical analysis using Jones' equation described 
above. Trends were evident but they were sometimes contrary to 
those predicted by Jones' analysis. The present experimental 
test program aims to determine which method, classical 
analysis or FEA, produces the best results, both in accurately 
predicting buckling strengths and in the relative ranking of 
different laminate types. 
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2.7. Southwell- type Analysis of Elastic Stability 

In 1932, R.V. Southwell, in an attempt to find "the 
practical value of a theory of elasticity, " proposed a 
procedure to analyze test data from an elastic column with 
initial curvature in an effort to determine the buckling 
characteristics of a perfectly straight column [6]. This 
procedure was an attempt to validate theories of elastic 
stability experimentally, realizing that "(1) the finite 
strength of actual materials, and (2) the unavoidable 
imperfections of workmanship" prevent the practical 
realization of elastic stability's concept of a "critical 
load." In the method, the central deflection of the column was 
plotted against a load-normalized deflection and the slope of 
the best straight line drawn through the test data was 
interpreted as the critical buckling load of the a perfectly 
straight column. Southwell foresaw that the basic concepts of 
this method should be widely applicable to structural 
stability problems wherever "inaccuracies in the specimen or 
in the experimental apparatus introduce displacements which 
increase continuously with the load." 

Recognizing its wide applicability to elastic stability 
problems, in 1955, Galletly and Reynolds developed an 
extension of Southwell's method that could be used to 
determine the elastic general instability pressure of ring- 
stiffened cylindrical shells subjected to external hydrostatic 
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pressure [7]. In this procedure the applied hydrostatic 
pressure was plotted against the circumferential bending 
strain that existed largely due to the initial out-of- 
roundness of the test specimens. 

L.H. Donnell concluded back in 1938 that Southwell's 
straight-line plot method was probably not applicable to 
curved structures because the deformed shape of curved 
structures was not a developable surface and the corresponding 
differential equations would be nonlinear [31]. Galletly 
and Reynolds showed that although the plot of pressure versus 
the reciprocal of the bending strain due to out of roundness 
was not linear, it was so slightly curved as the buckling 
pressure was approached that it could be closely approximated 
by a straight line. Another paper, by Horton and Cundari, 
demonstrated using previously published experimental data and 
theoretical analysis that the method was "of absolute 
generality, " and they concluded that "theories for perfect 
bodies . . . not in agreement with experimental values derived 
on the basis of Southwell type plots must be considered 
inadequate, incomplete or in error." [32] 

Galletly and Reynold's method was based on small- 
deflection theory where displacements and strains grow 
indefinitely large as the buckling pressure is approached. 
Because strains are easier to measure experimentally than 
displacements and intercepts are easier to objectively 
determine than asymptotic slopes, they plotted the reciprocals 
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of the strain quantities versus pressure rather than 
displacement versus pressure like Southwell had originally. 

An explanation of the method of applying Southwell's 
concepts that was used to analyze the strain data from this 
test program is presented graphically in the next three 
figures. Figure 2.5 shows some typical circumferential strain 
data plotted versus the applied hydrostatic pressure. The 
curved lines show the variation of the circumferential strain 
with pressure at three individual points on the circumference 
of a cylindrical specimen (0°, 45°, and 90° relative to each 
other). These curves become hyperbolic as the critical load is 
approached, but notice that they can be approximated by 
straight lines for a large portion of their length, as shown 
by the dashed lines. The slopes of these straight lines are 
different at each different location around the cylinder due 
to the varying amounts of local eccentricities. Galletly and 
Reynolds illustrated that the average of a large number of 
evenly distributed circumferential strain gage readings 
closely approximates the axisymmetric circumferential strain 
predicted by the Equation 13. This membrane strain can be 
considered as the elastic circumferential strain that would 
exist in a perfect structure whose radius is being uniformly 
compressed. The individual circumferential strains that 
actually occur on the cylinder are due to a combination of 
this membrane strain and the bending strain induced by the 
initial local eccentricities in the circular cylindrical shape 
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Figure 2.5. Typical plot of circumferential strain versus applied 
hydrostatic pressure, showing extraction of bending and nonlinear 
strains . 

of the structure. We therefore consider the average of the 
experimental strains as the membrane strain, G,e»brane, and the 
difference between this and the actual individual strain gage 
readings as the local bending strain, Gbending, as shown in 
Figure 2.5. 

A plot of the bending strain versus the applied pressure, 
p, is nonlinear and Gtenaing asymptotically approaches infinity 
at certain values of p, the first of which is the minimum 
elastic buckling pressure, Pcrit* The other values of p where 
^bending becomes infinite do not occur in reality and therefore, 
a plot of the reciprocal of the experimentally measured 
bending strain, 1/Gbending, versus the applied pressure will be 
curve that intercepts the p-axis at Pcit/ the minimum buckling 
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Figure 2.6. Plot of l/^bending versus pressure to determine a 
cylinder's critical buckling pressure. 

pressure of a perfect cylinder. 

Referring again to Figure 2.5, the gages located at the 
relative locations of 0“ and 45° behave in a manner that makes 
this method of determining Pcrit feasible. However, the fact 
that the strains from the gage located at 90° cross the 
membrane strain line prevents the utilization of data from 
this gage in this exact method. This is because where € equals 
the vslue of l/Gbending becomos infinite. The Southwell 
plot of all three gages where Gaverage = 6.e.brane is shown in 
Figure 2.6. 

One way to be able to use the readings from all "good" 
gages is to plot the reciprocal of only the nonlinear 
component of the local bending strain, instead of the bending 
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Figure 1 . 1 . Plot of 1/^non-linear versus pressure to determine a 
cylinder's critical buckling pressure. 



strain itself. In Figure 2.5, the strain at each gage is 
noticed to be very linear for the first 50 to 75 psi of 
pressure. The dashed straight lines are extensions of this 
linear portion of the strain readings. We can, therefore, 
break each bending strain response curve into a linear portion 
and a nonlinear portion we will call €ooniinear, such that; 



’^bending 



= Cp + € 



nonlinear 



(63) 



Where c is the constant slope of the dashed line. A plot of 
Gooniinear versus the pressure also, as expected, becomes infinite 
at the critical buckling pressure. The three gages of Figure 
2.5 were used in this manner to create the plot of l/G„oDiiDest 
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versus p in Figure 2.7. 

When one recognizes that the above deduction makes no 
assumptions or references to the material characteristics of 
the cylinder, the applicability to the difficult problem of 
laminated fiber reinforced composite cylinders is realized. 
These methods were used in an attempt to validate previous 
theoretical and finite element analysis of laminated composite 
cylinders, using strain data collected from the hydrostatic 
test of composite cylindrical specimens. Predicting the 
critical buckling pressure during a test before the onset of 
failure is also a possibility with this method. It is proposed 
that a method could be developed whereby strain data from a 
pressure hull loaded to only a percentage of its predicted 
buckling load could be used to predict how that cylinder will 
behave at pressures approaching the critical buckling 
pressure. Development of a method of this type would enable 
for hull certification and in-service monitoring to determine 
degradations from impact damage, cyclic loading, etc. 
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2.8. Buckling Mode Shape Analysis 



A method is proposed to determine the buckling mode shape 
of a circular cylinder using the basic concepts that R.V. 
Southwell presented for relating the radial displacement of a 
point on a axisymmetrically compressed cylindrical tube to its 
change in circumference [33]. For this analysis, the 
radial displacement is assumed to occur in a plane 
perpendicular to the longitudinal axis, and the circular 
cross-section is assumed to remain circular. The fractional 
contraction of the circumference results in a new radius (R- 
w), and can be expressed as an axisymmetric circumferential 
strain, 

2%{R-w)-2%R _ _w (64) 

y ■ 2tiR ■ R 

This can be converted into an expression for the radial 
displacement in terms of just the nominal radius R and the 
uniform circumferential strain 6y, 

w = -&yR (65) 

The actual circumferential strain at any point on a cylinder 
being compressed nonuni formly can be broken into a membrane 
stain that would define the uniform contraction in the radius 
and a bending strain that would be the local contraction or 
extension relative to this uniform change. Therefore, by 
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looking at just the affect of the bending strain at several 
points around the circumference, an idea of the number of 
circumferential buckling waves can be deduced. 

This method is applied to each of the specimens using the 
eight circumferential strain gage readings at the last 
pressure where all gages were recorded prior to buckling 
failure. To better illustrate the deformed shaped the radial 
deflections were magnified from 200 to 500 times normal size. 

These mode shape plots are contained in Appendix C. The 
odd-numbered figures contain the variations of bending strain 
with pressure at the noted locations around the circumference 
of each specimen and the even-numbered figures contain polar 
plots of the deduced mode shapes with the magnified radial 
deflections . 

Application of the above analysis to these bending 
strains and magnifying the radial displacements, even though 
it does not give absolute deflections, does make the buckling 
mode shape clearly apparent. This experimentally derived mode 
shape can then be used directly in the classical and finite 
element buckling analyses, or as a verification of the 
analyses' circumferential buckling wave assumptions. 
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3. EXPERIMENTAL PROCEDURE 



The main contribution of this thesis is the results 
obtained from the hydrostatic pressure testing to failure of 
cylindrical test specimens of various laminate configurations. 
The purpose of this testing was to experimentally verify the 
results obtained from classical and finite element buckling 
analyses. In addition to checking how precisely the analysis 
methods predict the buckling failure pressures of these 
relatively simple laminates, these methods must be verified to 
at least rank the buckling pressures of the different laminate 
types correctly. The correct ranking is essential to selecting 
one of these analysis methods as a basis for studies into the 
optimization of ply orientations for buckling strength. 

The test matrix originally called for the testing of at 
least three 12-ply composite cylinders from each of the four 
laminate configurations. In all, thirteen cylinders were 
tested but because of fabrication induced defects only eight 
specimens were considered "good." The additional 5 defective 
specimens were tested to failure, but their strain data and 
comparative value were rejected. These cylinders had 
noticeable flaws and would normally have been discarded prior 
to testing, but in order to refine the test set-up and also in 
an effort to gain an understanding of how these defects affect 
the failure mechanisms, all were tested, but only the strain 
data from the good specimens was used in the analysis and 
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presented in the results. 

The strain gage readings taken throughout the hydrostatic 
tests and the critical buckling loads measured comprise the 
bulk of the experimental data. The results were tabulated and 
plotted in various formats, and then analyzed to determine 
patterns or conclusions about the strength of composite 
cylinders subjected to external hydrostatic pressure. The 
experimental results were also compared with the analytical 
results to check predictions about laminated cylinder buckling 
strengths and assumptions of buckling mode shapes. 

Very little research has been conducted to experimentally 
determine the behavior of thin-walled composite cylinders of 
various lay-ups subjected to hydrostatic external pressure to 
failure, and to compare the results with theoretical 
predictions. Currently there is only a meager amount of test 
data on the buckling of thin-walled composite shells under 
hydrostatic pressure. The data collected and analyzed here is 
hoped to give a better understanding of the ability to use 
composites for submersible hull structures. Ultimately, 
insight is sought into improved ways to select material 
combinations and laminate configurations that will enable the 
benefits of laminated composite material technology to be 
maximized in future design and construction of submersible 
pressure hulls. 
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3.1. Cylindrical Test Specimens 



All of the tested specimens were fabricated from 12-ply 
laminated circular cylinders with nominal diameters of 5.52 
inches, skin thicknesses of 0.072, and lengths of 18 inches. 

All the cylinders were constructed of individual plies of 
T300/5208 graphite/epoxy pre- impregnated unidirectional ply 
(prepreg uniply). Table 3.1 contains the nominal properties of 
T300/5208 in this format and the typical properties of various 
other types of graphite fiber / polymer matrix unidirectional 
composites, for comparison. These properties were obtained 
from the Composite Materials Research Group at the University 
of Wyoming and are for nominal 60% fiber volume composites. 
Digital micrographs of polished specimens cut from the excess 
ends of the laminated tubes tested here indicated an average 
fiber volume of 50%. 
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Composite 


AS4/ 


GY70/ 


T300/ 


Property 


3501-6 


930 


5208 


E„ 


(Msi ) 


20.0 


47.1 


19.0 


E22 


(Msi ) 


1.3 


0.9 


1.5 


Gi2 


(Msi ) 


1.0 


0.6 


1.0 


V 12 


0.30 


0.26 


0.28 


ail 


(10‘V°C) 


-0.04 


00 

• 

0 

1 


-0.04 


a 2 2 


(10‘V°C) 


18.0 


18.0 


18.0 


a.t 


(ksi) 


210 


110 


218 


a2t 


(ksi) 


7.5 


3.8 


5.8 


Ti2 


(ksi ) 


19.0 


3.9 


13.9 


a 1 c 


(ksi) 


210 


102 


218 


O20 


(ksi ) 


29.9 


5.0 


35.7 



Table 3.1. Typical Material Properties. 



A. Fabrication Method 

There are several different viable fabrication methods 
that could have been used to produce the cylindrical specimens 
for the hydrostatic testing. The primary methods used in the 
construction of tubes and cylindrical vessels out of laminated 
composite materials are hand layup with vacuum bag curing, 
filament winding, and tape winding. Hand layup with vacuum bag 
curing was not considered to be adequate because it would 
almost be impossible to control wrinkles and geometric 
deviations from perfect circularity. As discussed earlier, 
slight imperfections have a significant effect on the buckling 
failure characteristics of cylindrical shells in compression. 
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Alum Upper External Tool 





Laminated Cylindrical Specimen 



Figure 3.1. The elastomeric tooling fixture used to manufacture the 
cylindrical test specimens. 



These irregularities could be largely eliminated and at least 
controlled by constructing the specimens with filament winding 
technique. This fabrication method produces highly circular 
and wrinkle-free cylinders, but distribution of the fibers 
could not have been made as uniformly as required. Tape 
winding, which combines the best attributes of hand layup and 
filament winding was a technically feasible option, but the 
cost of a tape winder was prohibitive to project budget 
constraints. For these reasons, a novel method of constructing 
the cylindrical specimens with high circularity and minimal 
laminate defects had to be incorporated. An elastomeric 
tooling method that involves the use of an expanding internal 
mandrel and a highly circular external mandrel was chosen to 
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fabricate the test specimens. 

This method calls for an internal mandrel constructed of 
a high coefficient of thermal expansion (CTE) material which 
expands in radius during the increased temperature curing 
cycle. This mandrel is wrapped with thin sheets of elastomer 
to give a consistent, smooth surface to the final cylinder. 
Sheets of prepreg uniply are then hand laid up onto the 
mandrel. The sheets of uniply must be carefully measured, cut 
and rolled onto the cylinder by one person while the another 
person makes sure that the lamina lays down smoothly and at 
the desired angle. After all plies are applied to the internal 
mandrel, two highly accurately machined semi-circular metal 
external tools are clamped on. The pressure exerted during the 
curing process by the expanding high CTE internal mandrel 
forces the multiple layers of uniply tape up against the 
stiff, low CTE external tool. 

This method proved to consistently create laminated 
cylindrical specimens with extremely close circular 
tolerances. The details of the elastomeric tooling fixture 
used to manufacture the specimens for this research is 
diagrammed in Figure 3.1. It has an internal mandrel made of a 
syntactic foam core (CTE = 28 x 10'® in/in°F) that is wrapped 
with several layers of silicone rubber sheet (CTE = 162 x 10'® 
in/in°F), and an external tool fabricated from two 1/2" thick, 
24" long aluminum pipe halves (CTE = 14 x 10'* in/in°F) which 
were bored to a .005" circular tolerance. 
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B. Laminates Tested 



The cylinders each had twelve plies and were constructed 
out of one of the four following different layups: 

[904/02]s - 2:1 specially orthotropic 
[QOa/Osls - 1:1 specially orthotropic 
[902/(±30)2 ]s - quasi-isotropic 

[902/0]«t - interleaved orthotropic 
The numbers indicate the orientation of the individual plies, 
the numerical subscripts indicate the number of plies of that 
orientation, the alphabetic subscript S indicates that the 
laminate is symmetrical about the mid-surface, and the 4T 
indicates that the unsymmetric laminate is constructed of 4 
identical plies. Each of these laminates are diagrammed in 
Figures 2.2, 2.3, 2.4, and 2.5. 
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Figure 3.2. 
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Figure 3.3. 
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Figure 3.4. 



Figure 3.5. 
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The 1:1 and 2:1 specially orthotropic laminates do not 
exhibit any stiffness coupling and, therefore, are the easiest 
of the anisotropic laminates to analyze. This also means that 
the material does not warp or twist during post-cure cooling. 
This is not as much of a problem with cylinders as flat 
plates, but was a concern nonetheless. 

The unsymmetric, interleaved orthotropic layup, which is 
commonly used in industry, has somewhat better interlaminar 
properties than the specially orthotropic laminates because 
the perpendicular ply interfaces are distributed through the 
thickness. This causes a reduction of the interlaminar 
stresses that cause delamination, a major failure mechanism in 
laminated composites. 

The quasi-isotropic layup was used to compare the 
composite cylinder failure characteristics to cylinders made 
from isotropic materials. Laminates made with this type system 
of ply orientations result in directional stiffness 
characteristics very close to isotropic. 
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3.2. Hydrostatic Pressure Tests 



A. Test Set-up and Instrumentation 

The hydrostatic testing was conducted in an Atlas 600 psi 
air pressure chamber located in the Mechanical-Thermal Lab at 
C.S. Draper Laboratory in Cambridge, MA. Pressure was 
provided to the chamber via a pressure control system that was 
accurate to within ±2 psi. The cylinders were potted at both 
ends in specially designed endcaps. The endcaps, shown in 
figure 2 . 6 , were designed with a slot for potting that has an 
angled outside edge so that the edges remained circular, and 
also so that only the outside edge of the tube would actually 
be in contact with the endcap. Armstrong C7 epoxy was used to 
pot the tubes in the endcaps. This epoxy was compliant enough 
to allow the ends of the cylinders to rotate when loaded, 
simulating simply supported edge conditions. 

Each specimen was instrumented along the centerline with 
10 Micro Measurements ^5 inch, 120 Q strain gages, model number 
EA-06-500BH-120 . Eight of the gages were evenly distributed at 
45° separations around the centerline of the specimen, and 
oriented to read circumferential or hoop strain. The other two 
gages were positioned on the centerline opposite each other at 
a 180° separation, and oriented to read axial strain. The gage 
leads were draw out of the pressure chamber, via breakout 
boxes, to a ten channel BLH Model 1225 switching and balancing 
unit, and then to a BLH Model 1200 digital strain indicator. 
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The test set-up is diagreunmed in Figure 3.7 and a photograph 
of the facilities is provided in Appendix E. 

B. Test Procedure 

The testing of each cylindrical specimen was conducted by 
the following procedure: 

o The specimen was placed in the chamber, the strain gage 
leads were connected to the breakout box, the gages tested for 
integrity, and a copper bleed-off line was connected to one of 
the endcaps. 

o The pressure chamber was closed and sealed, and each of 
gages was manually zeroed with the switching and balancing 
unit . 

o The pressure in the chamber was increased to 25 psi and 
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Pressure Chamber 





Figure 3.7. Diagram of experimental set-up, showing pressure chamber 
and instrumentation. 

maintained long enough for strain readings to stabilize and to 
be manually recorded for each individual gage, 
o The pressure was increased to 50, 75, 100, 110, and so on at 
10 or 5 psi increments until failure, taking readings at each 
step. 
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3.3. Ring Tests 



The external tools that were used in fabrication of the 
cylinders were 24 inches long but the specimens were only to 
be 18 inches long. This left 3 inches of excess on each end of 
the cylinder. Several 1/2 inch rings were cut from the excess 
end of each cylinder, and these were used to help determine 
some of the in situ material properties of the specimens. 

Two strain gages were applied at the same angular 
location of each ring, one on the inside and the other on the 
outside surface, 90° from the point of load application. The 
rings were then loaded in two point compression on the 
Instron® 4505 located in the Composite Materials Laboratory of 
the Material Science and Engineering Department at the 
Massachusetts Institute of Technology, Cambridge, MA. 
Measurements of strain and displacement versus load were 
recorded throughout each of the tests. The data acquisition 
system consisted of a four channel Micro Measurements strain 
gage conditioner unit, model 1220A, and a National Instruments 
Lab-PC integrated data acquisition card. 

These ring tests gave some additional mechanical property 
data for the various laminate configurations in the "as 
fabricated" condition. 



74 



4. RESULTS and DISCUSSION 



4.1. Ring Tests 

Throughout the two-point, symmetric compression test of 
the each 1/2 inch ring, load, strain and displacement data was 
collected. Plots were made of the load versus displacement and 
the load versus strain data. The slopes of the load versus 
displacement curves were used with the relationship of 
equation 28 to calculate a "w test" for each of the ring 
specimens. Additionally, the slopes of the strain versus 
displacement curves were used with the relationship of 
equation 31 to calculate an "e test" D 22 . These experimentally 
determined D 22 ' s were compared with the D 22 from CLT contained 
in the ABD matrices of Appendix E, which were developed using 
the material properties for T300/5208 given in Table 3.1. This 
comparison is summarized in Table 4.2 below. The experimental 
"w-test" and "€-test" values are seen to agree well with each 
other, and are in reasonable agreement with the values 
calculated using CLT. 

Because the radius, thickness and width of the ring 
specimen has a significant impact on equation 31, these 
dimensions were accurately measured with a micrometer at eight 
locations around the ring and the average of these values used 
in calculations. These averaged values for R, t and b are 
presented in Table 4.1. 
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Specimen 

No. 


R 

[in. ] 


t 

[in.] 


b 

[in. ] 


6 


2.759 


0.072 


0.523 


9 


2.761 


0.070 


0.503 


12 


2.757 


0.074 


0.516 


13 


2.759 


0.072 


0.505 


17 


2.758 


0.077 


0.510 


18 


2.757 


0.076 


0.506 


23 


2.760 


0.075 


0.507 



Table 4.1. Dimensions of Ring Specimens. 



Specimen 

No. 


P/w 


P/e 


D 22 

w-test 


D 22 

G-test 


D 22 

CLT 


6 


147.1 


12650 


440.1 


436.6 


526.2 


9 


116.3 


10100 


361.6 


352.1 


526.2 


12 


128.3 


10520 


387.5 


377.4 


574.4 


1 


134.9 


11180 


417.5 


399.6 


574.4 


17 


126.0 


9278 


385.9 


351.1 


445.7 


18 


147.0 


11400 


453.3 


428.9 


445.7 


23 


93.4 


8623 


288.5 


319.9 


407.0 



Table 4.2. Predicted and Measured D 22 
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4.2. Cylinder Tests 



A total of 13 full-length cylindrical specimens were 
tested to failure, but several of these had noticeable defects 
and thus only served to refine the test procedure and gain 
more experience with the graphite/epoxy pressure vessels with 
large defects. Table 4.3 contains a breakdown of the specimens 
tested and their critical buckling pressures. This table also 
gives comments on each "bad" specimen, explaining the details 
of the problem, such as any noticeable flaws. In all, eight 
specimens have been tested that were classified as "good". 

Only the buckling pressures and strain data from these eight 
specimens were analyzed for this investigation. 



# 


Type 


P crit 


Notes and Comments 


Class 


5 


1:1 


140 


Tested with conical steel endcaps. 


bad 


6 


1:1 


183 




good 


7 


1:1 


95 


Creased when clamped into tool. 


bad 


8 


1:1 


160 


Creased when clamped into tool. 


bad 


9 


1: 1 


200 




good 


10 


2:1 


167 


Coated in epoxy due to leakage. 


bad 


12 


2:1 


181 




good 


13 


2:1 


173 




good 


15 


Unsym 


67 


[02/±604/02]t. Used steel endcaps. 


bad 


16 


Quasi 


165 




good 


17 


Quasi 


183 


Loaded 50 psi, unloaded, reloaded. 


good 


18 


Quasi 


175 




good 


23 


Unsym 


152 




good 



Table 4.3. Summary of all tested specimens. 
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The original endcaps were designed with large cones to 
seat the cylinder edges onto and were constructed out of 
steel. Reproducible end conditions were extremely difficult to 
create using these steel endcaps. Specimens 5 and 15 were the 
first two cylinders tested and the conical endcaps were used. 
The cones on the endcaps seemed to wedge into the cylinders 
under pressure creating erroneous results. These problems 
inspired the development of a new endcap design. The details 
of the potting of the cylinders in these endcaps are shown in 
Figure 3.6. Because of the change in endcaps the results from 
the first two cylinders tested were discarded. 

Specimen 10 was a very "dry" laminate and had excessive 
leakage through the plies themselves. This specimen was coated 
in C-7 epoxy to prevent the pressure from equalizing inside 
the cylinder and enable testing to be completed. Because of 
this coating, the results were not considered accurate. 

During the manufacture of specimens 7 and 8 the pre-cured 
laminate that was rolled onto the internal mandrel was too 
large for the aluminum external tools to be clamped around. 
These specimens were put into a freezer to induce shrinking 
and then forcibly clamped into the tools. Both these specimens 
had creases of the order of the laminate thickness from the 
pinching of the semicircular aluminum tools. Failure of these 
specimens under hydrostatic loading was premature and both 
fractured cleanly down the creases as shown in the photographs 
in Appendix F. For this reason these specimens were not used 
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for comparative purposes. 

Table 4.4 summarizes the critical buckling pressures that 
were obtained from Jones' equation, ABAQUS modelling, 
experimental testing, and Southwel 1-type analysis for the 
"good" specimens. 

Although the test matrix was not as large as originally 
planned, "good" data was obtained for two 1:1 specially 
orthotropic, two 2:1 specially orthotropic, three quasi- 
isotropic, and one interleaved ( unsymmetric ) orthotropic 
specimen. Because of the close similarity between the data 
from each of the "good" specimens from the same layup, these 
results can be confidently used to evaluate the analysis 
methods . 



Layup 


Cyl inder 


Jones ' 


ABAQUS 


Test 


Southwell 


[903/03]s 


6 


214 


194 


183 


200 


[903/03]s 


9 


214 


194 


200 


200 


[904/02]s 


12 


196 


204 


181 


188 


[904/02]s 


13 


196 


204 


173 


177 


[902/(±30)2]s 


16 


185 


172.5 


165 


167 


[902/(±30)2]s 


17 


185 


172.5 


183 


183 


[902/(±30)2]s 


18 


185 


172.5 


175 


179 


[902/0]4T 


23 


164 


197 


152 


155 



Table 4.4. Predicted and Measured Critical Buckling Pressures. 



A. Critical Buckling Pressures 

As stated earlier, a major objective of this work is to 
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verify which of the analysis methods not only best predicts 
the buckling characteristics, but consistently correctly ranks 
the stability of each of the laminates used. This latter 
objective is crucial as these methods are proposed for further 
studies into optimizing the orientations of the plies in a 
multi-layered composite cylinder to achieve the highest 
buckling strength. 

Figures 4.1, 4.2, 4.3 and 4.4, contain bar graphs which 
each present the ranking of the four different laminate types. 
Figure 4.1 gives the values obtained using Jones' buckling 
criterion, equation 56, which was derived from Donnell's 
stability equations and CLT. Figure 4.2 are the values given 
by the FEA program ABAQUS. Figure 4.3 presents the average 
critical pressure for each laminate type predicted by the 
Southwell -type analysis, taken from the plots in Appendix B. 
Figure 4.4 gives the average experimental buckling pressures 
for each laminate type as obtained from the hydrostatic tests 
of the "good" specimens. 

These figures show that the ABAQUS finite element 
analysis program does not rank the laminates correctly. This 
is a major obstacle to using ABAQUS to optimize the selection 
of individual ply orientations. The ABAQUS analysis predicts 
that the orthotropic laminates with two hoop plies for every 
axial ply will have the greatest buckling strength. It 
predicts that the 2:1 symmetric specially orthotropic laminate 
is strongest and that the interleaved, unsymmetric orthotropic 
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laminate will have a higher critical buckling pressure than 
either the 1:1 specially orthotropic of the quasi-isotropic 
laminates. This was not collaborated by the testing. 

The Jones' equation analysis predicted that the 1:1 
specially orthotropic laminate would give the cylinder the 
greatest buckling strength, followed by the 2:1 specially 
orthotropic, the quasi-isotropic, and then the interleaved 
orthotropic laminates, respectively. This ranking was exactly 
the same order as the critical pressures measured in 
hydrostatic tests and predicted by the subsequent Southwell- 
type analysis. This is a significant finding. Applying 
optimization techniques to the classical analysis methods of 
Jones will be much less expensive than to apply similar 
techniques to the finite element methods of the ABAQUS 
program. That these classical methods have been shown to 
correctly rank at least the popular laminate types tested 
here, means that designers can apply them to optimization 
studies with a substantial amount of confidence. 
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Figure 4.1. Laminate ranking using Jones' Equation. 
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Figure 4.2. Laminate ranking using ABAQUS. 
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Southwell Analysis Predictions 
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Figure 4.3. Laminate ranking using Southwell- 
analysis . 
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Figure 4»4. Laminate ranking using the results 
hydrostatic testing. 
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Next, the actual critical pressures from the hydrostatic 
tests are compared to the buckling pressure predictions to 
verify the accuracy of the various analysis methods. Figures 
4.5, 4.6, 4.7 and 4.8 contain bar graphs which each present 
one laminate type and show a quantitative comparison of the 
different buckling pressures predicted by the three analysis 
methods and obtained from hydrostatic testing. 

The ABAQUS program was very accurate in its prediction 
for the cylinders constructed out of the 1:1 specially 
orthotropic and the quasi-isotropic laminates, but the 
magnitude of the overestimation of the buckling strength for 
the 2:1 specially orthotropic and interleaved, unsymmetric 
orthotropic laminate cylinders prevents this finite element 
method from being used with any confidence to predict the 
critical pressures of angle ply laminated cylinders. 

Jones' equation gives predictions of the buckling 
strength larger than those measured experimentally, but within 
reasonable deviation. It is significant to note the precision 
and accuracy with which the Jones' equation and the average 
Southwell -type analysis exhibit, relative to each other. Table 
4.5 contains the results of these two analyses which both 
consider the strength of a "perfect" geometry circular 
cylinder. It should be noted that the Jones' predictions were 
for 60% fiber volume laminates, and recently taken digital 
micrographs of the specimens indicate an average 50% fiber 
volume. Recalculations using the updated laminate properties 
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show that the values obtained from the two methods are 
actually in even better agreement. 



Laminate 


Jones ' 


Southwel 1 


Percent 


Type 


Equation 


Plot 


Difference 


1:1 Specially Orthotropic 


214 


200 


+ 7.0 % 


2:1 Specially Orthotropic 


196 


182.5 


+ 7.4 % 


Quasi- isotropic 


185 


176 


+ 5.1 % 


Interleaved Orthotropic 


164 


155 


+ 5.8 % 



Table 4.5. Comparison of Jones' Equation and Southwell Plot. 
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1:1 Specially Orthotropic 




Figure 4.5. Critical buckling pressures for 1:1 Specially 
Orthotropic laminates. 



2:1 Specially Orthotropic 




Jones’ Eqn 



ABACUS 



Southwefl 



Test Avg 



Figure 4.6. Critical buckling pressures for 2:1 Specially 
Orthotropic laminates. 
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Quasi-isotropic 




Figure 4.7. Critical buckling pressures for Quasi- 
isotropic laminates. 



Interieaved Orthotropic 




Figure 4.8. Critical buckling pressures for Interleaved 
Orthotropic laminates. 
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B. Axial Strain Response 

It was hoped that the two axially mounted strain gages 
could be used to investigate the beam-column effect of the 
long cylindrical specimens. The plots of the axial strain 
responses for all of the good specimens are contained in the 
even-numbered figures of Appendix A. The axially strains for 
the 1:1 specially orthotropic specimens seemed to be quite 
linear with little or no increase in curvature as the critical 
pressure was approached. The plots of the axial strain 
response for the 2:1 specially orthotropic, quasi-isotropic, 
and interleaved (unsymmetric) orthotropic specimens showed the 
more traditional hyperbolic shape at increased pressures. 

C. Circumferential Strain Response 

Plots of the circumferential strain responses from the 
strain gages located on each of the specimens are contained in 
the odd-numbered figures of Appendix A. All the cylinders had 
eight circumferentially-mounted strain gages evenly 
distributed around the centerline at 45° separations. Note 
that when one of the plots has less than eight circumferential 
strain curves, it means that one or more of the gages failed 
during the test. 

The hyperbolic curvature of the circumferential strain 
response measured by each of the gages, as predicted by 
theory, is clearly evident in all cases. The individual 
circumferential strain responses to increasing hydrostatic 
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pressure can be broken into a membrane and a bending strain 
response. The bending strain results from the shape change of 
the circular cylinders' cross-section as the buckling 
pressure is approached. 

This excellent circumferential strain data is a 
significant contribution of this research, and proved well- 
worth the time-intensive labor it required to apply the gages, 
and then to collect and analyze the extensive amount of data 
that was yielded. Because the eight circumferential strain 
gages were evenly distributed around each cylinder's 
centerline, they provided the necessary data to construct 
detailed Southwel 1-type plots and to perform an analysis to 
deduce the buckling mode shape for each cylinder. 

D. Average Circumferential Strain 

By averaging the readings of the eight circumferential 
gages at each pressure we get a close approximation of the 
axisymmetric circumferential membrane strain of each of the 
cylinders as they were compressed. This is the strain response 
that would exist if the cyl inders* cross-sect ion would remain 
circular throughout the loading. 
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Figure 4.9. Axisynunetric circumferential membrane strain 
response for the 1:1 specially orthotropic cylinders. 




Figure 4.10. Ax i symmetric circumferential membrane strain 
response of the 2:1 specially orthotropic cylinders. 
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Figure 4.11. Axisymmetric circumferential membrane strain 
response for the quasi-isotropic cylinders. 




Figure 4.12. Axisymmetric circumferential membrane strain 
response for the interleaved orthotropic cylinder. 
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To examine how accurately the averages of the eight 
circumferential strain gage readings approximate the membrane 
strain, the experimental average of the circumferential strain 
gage readings is compared to the theoretical axisymmetric 
circumferential membrane strain response predicted by equation 
9 for hydrostatic loading. The curves for each specimen and 
the theoretical curve for each laminate type, are contained in 
Figures 4.9, 4.10, 4.11 and 4.12. 

These plots show that the average of the eight 
circumferential strain gage readings from specimens made from 
the same laminate, are almost exactly the same. These average 
strains also compare extremely well with the axisymmetric, 
circumferential membrane strain calculated using the material 
properties for T300/5208 with 60% fiber volume and equation 9, 
which is based on CLT. The slight discrepancy could be the 
result of the fact that the laminates actually had only a 50% 
fiber volume, not 60% as assumed in the material properties. 
Additionally, equation 9 does not consider the bending 
response or the extension/shear coupling in the laminate. 

E. Modified Southwell Analysis 

Using the method described in the section on Southwell 
analysis, where the reciprocal of the non-linear strain has 
plotted versus pressure, Southwell -type plots were constructed 
for each good specimen. The odd-numbered figures of Appendix B 
show the Southwell plots for each of the specimen over the 
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same large range of pressure. The even-numbered figures of 
Appendix B show the details of the Southwell method by 
displaying only a narrow range of concern for each specimen. 
The dotted lines on these plots are the extensions of the 
reciprocals of the non-linear strain curves. Close exeimination 
of these plots show the extraordinary small amount of scatter 
that was involved in the prediction of critical buckling 
pressures. The values of these predicted buckling pressures 
for the specimens had they been of perfect geometry is given 
in the tables and bar charts earlier in this chapter. As noted 
before, the predictions of Jones' equation were consistently 
only between 5 to 7.5 percent higher than those obtained from 
the Southwell plots. 

F. Mode Shape Deduction from Bending Strains 

Using the method outlined in the section on Mode Shape 
Analysis, the bending strains of each specimen were utilized 
to deduce the buckling mode shape of each cylinder. The plots 
of the bending strains and the subsequent variations in shape 
are contained in Appendix C. The buckling mode shapes are 
drawn at the last pressure where all the strain gages were 
read and recorded. 
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Layup 


Specimen 


Mode Shape Deduction 






Test 


CLT 


FEA 


[903/03]s 


6 


3 


3 


2 


[903/03]s 


6 


3 


3 


2 


[90,/03]s 


12 


3 


3 


2 


[90,/03]s 


13 


3 


3 


2 


[903/(±30)3]s 


16 


3 


3 


3 


[903/(±30)3]s 


17 


3 


3 


3 


[903/(±30)3]s 


18 


3 


3 


3 


[903/0]4t 


23 


3 


3 


3 



Table 4.6. Predicted and Experimentally Deduced Mode Shapes. 

It is significant that each and every specimen, 
regardless of laminate type, appeared to buckle in a 
distinctive n = 3, three lobe mode shape. This was predicted 
by Jones' equation for all cases, whereas the FEA predicted 
n=2 to be the lowest energy buckling mode shape for both the 
1:1 and 2:1 specially orthotropic cases. 

G. Failure Progression 

The 5208 epoxy matrix used in the laminates tested in 
this investigation is a very brittle resin, having a low 
toughness and high modulus (1.4% failure strain). Brittle 
resins tend to be unable to resist the propagation of 
delamination, exhibiting quite sudden and uncontrollable 
failures . 

The failure of these T300/5208 laminated cylinders was 
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exceptionally catastrophic, appearing to happen in 
microseconds. Attempts were made to film the collapse of the 
cylinders with a high-speed (250 frames per second) video 
camera, but even when slowed to one frame per second for 
viewing the failure sequence was too fast to capture on film. 

All the good specimens shattered spectacularly, leaving 
only very small pieces of intact laminate, large quantities of 
individual broken fibers, and basically only a dust residue of 
the resin. Both the specimens that were creased during 
fabrication (cylinders 7 and 8) fractured cleanly along the 
crease, leaving large sections of the cylinder relatively 
intact . 
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5. SUMMARY AND CONCLUSIONS 



The hydrostatic tests conducted for this thesis have been 
successful in achieving the original objectives. 

The classical and finite element buckling analysis 
methods were examined for accuracy and correctness in 
determining the critical buckling pressure of unstiffened 
laminated composite cylinders of various layups. From this 
examination, Jones' equation for the buckling of generally 
orthotropic cylinders has proven accurate in predicting the 
buckling loads and correct in the ranking of different 
laminate types. The validity of this equation which has been 
experimentally verified, and its relative simplicity make it 
well-suited for use in optimization programs that could be run 
economically on personal computers. 

A systematic test program was developed whereby 
cylindrical specimens can be hydrostatically tested to 
buckling failure under repeatable conditions. The 
instrumentation used yielded valuable strain data that was 
used in Southwell-type analysis and also in determining the 
buckling mode shape of the specimens. 

Comparison of the Southwell plots and the Jones' equation 
predictions, showed the consistency with which this type 
analysis of circumferential strain data from the test of 
geometrically imperfect specimen can be used to predict the 
buckling strength of perfect cylinders. The applicability. 
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which had previously been shown for isotropic ring-stiffened 
cylinders, has proven to be applicable to unstiffened 
cylinders constructed of laminated anisotropic materials as 
wel 1 . 

Through careful documentation of the fabrication process 
and through the utilization of a large diversity of material 
characterization tests, insight into the role fabrication 
induced imperfections, residual and interlaminar stresses play 
in the degradation of pressure hull performance has been 
gained. In addition, ring compression tests have shown that 
the bending stiffnesses predicted by classical laminate theory 
may for symmetric laminates be exaggerated. 
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6. RECOMMENDATIONS 



This testing and analysis has shown the large increase in 
knowledge that can be gained by the extensive strain gage 
instrumentation. It is recommended that in future testing the 
pressure be steadily increased and the strains be recorded 
continuously on a computer data acquisition system. In these 
tests, the pump on the chamber would kick on and off with 
quite violent jolts attempting to maintain a steady pressure. 
This disrupted the stabilization of the strain gages and 
caused delays before readings could be recorded. The effect of 
these shocks and delays on the tubes, which is unquant if iable, 
would be eliminated with continuously increasing pressure 
loading . 

It is believed that by modifying the Southwell and mode 
shape analysis procedures presented here that a practical 
method could be developed whereby the compressive buckling 
strength of individual pressure hulls could be qualified. A 
properly instrumented pressure hull, tested to only a 
percentage of its predicted collapse load, would exhibit 
strains and displacements that may enable the non-destructive 
verification of predicted collapse loads. 

This test series should be extended to angle-ply 
laminates to check Jones' equation’s accuracy for non- 
orthotropic ply orientations. Notwithstanding the lack of 
angle-ply testing conducted here, it is recommended that 
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Jones' equation be utilized in laminate optimization 
investigations in lieu of finite element methods, which have 
shown to improperly rank the buckling strength of even the 
simple laminates tested here. 

Graphite/epoxy laminated thin cylindrical shells such as 
those tested here, have shown to have potential in the design 
of small submersibles that are only required to operate at 
shallow depths. It is recommended that in order to avoid 
catastrophic failure of the hull structure, and subsequent 
complete loss of the vehicle, that utilization of tougher 
resins be investigated. 
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APPENDIX A 



CIRCUMFERENTIAL AND AXIAL STRAINS 
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Figure A, 1 . Circumferential Strains around the centerline 
of cylinder 6. 




Figure A. 2. Axial Strains at the center of cylinder 6, 



103 




Figure A. 3. Circumferential Strains around the centerline 
of cylinder 9. 




Figure A. 4. Axial Strains at the center of cylinder 9. 
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Figure A. 5. Circumferential Strains around the centerline 
of cylinder 12. 




Figure A. 6. Axial Strains at the center of cylinder 12. 
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Figure A. 7. Circumferential Strains around the centerline 
of cylinder 13. 




Figure A. 8. Axial Strains at the center of cylinder 13. 
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Figure A. 9. Circumferential Strains around the centerline 
of cylinder 16. 




Figure A. 10. Axial Strains at the centerline of cylinder 
16. 
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Figure A. 11, Circumferential Strains around the centerline 
of cylinder 17. 




Figure A . 12 . Axial Strains at the center of cylinder 17. 
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Figure A. 13. Circumferential Strains around the centerline 
of cylinder 18. 




Figure A. 14 . Axial Strains at the centerline of cylinder 
18. 
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Figure A. 15. Circumferential Strains around the centerline 
of cylinder 23. 




Figure A. 16 . Axial Strains at the center of cylinder 23. 
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APPENDIX B 



SOUTHWELL-TYPE ANALYSIS OF STRAINS 



III 




Pressure (psi) 

Figure B . 1 . Southwell Plot for Cylinder 6. 




Pressure (psi) 

Figure B.2. Detail of Southwell Analysis for Cylinder 6. 



112 




Figure B.3. Southwell Plot for Cylinder 9. 




Pressure (psi) 

Figure B.4. Detail of Southwell Analysis for Cylinder 9. 
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Pressure (psi) 

Figure B.5. Southwell Plot for Cylinder 12. 




Pressure (psi) 

Figure B.6. Detail of Southwell Analysis for Cylinder 12. 
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Figure B.7. Southwell Plot for Cylinder 13. 




Figure B.8. Detail of Southwell Analysis for Cylinder 13. 
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Figure B.9. Southwell Plot for Cylinder 16. 
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Figure B.ll. Southwell Plot for Cylinder 17. 




Figure B.12. Detail of Southwell Analysis for Cylinder 17. 
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Figure B.13. Southwell Plot for Cylinder 18. 




Figure B. 14 . Detail of Southwell Analysis for Cylinder 18. 



118 




Figure B. 15 . Southwell Plot for Cylinder 23. 




Figure B.16. Detail of Southwell Analysis for Cylinder 23. 
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APPENDIX C 



BENDING STRAINS AND MODE SHAPE PLOTS 
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Figure C . 1 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 6. 



Cylinder No. 6 

OOj/Ojls 




♦ Radial deflections shown ot 1 80 psi, magnified 400 times. 

Figure C . 2 . . Shape variation due to bending strain for 
Cylinder 6. 



121 




Figure C.3. ♦ Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 9. 



Cylinder No. 9 

[903/03]s 




♦ Radial deflections shawn at 1 90 psi, magnified 400 times. 

Figure C. 4 . . Shape variation due to bending strain for 
Cylinder 9. 
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Figure C. 5 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 12. 



Cylinder No. 1 2 

ooyo^ls 




♦ Radial deflections shown at 1 80 psi, magnified 200 times. 

Figure C. 6 . . Shape variation due to bending strain for 
Cylinder 12. 
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Figure C.7. . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 13. 



Cylinder No. 1 3 

[goyo^is 




♦ Radial deflections shown at 1 70 psi, magnified 200 times. 

Figure C. 8 . . Shape variation due to bending strain for 
Cylinder 13. 



124 




Figure C . 9 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 16. 



Cylinder No. 1 6 

[902/(t30)2]s 




*Radial deflections shown at 1 60 psi, magnified 300 times. 

Figure C. 10. . Shape variation due to bending strain for 
Cylinder 16. 
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■'■i 







I 




Figure C. 11 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 17. 



Cylinder No. 1 7 

[902/(t30)2Js 




♦ Radial deflections shown at 1 80 psi, magnified 200 times. 

Figure C . 12 . . Shape variation due to bending strain for 
Cylinder 17. 
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Figure C. 13 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 18. 



Cylinder No. 1 8 

[902/(t30)2]s 




♦ Radial deflections shown at 1 70 psi, magnified 200 times. 

Figure C. 14 . . Shape variation due to bending strain for 
Cylinder 18. 
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Figure C. 15 . . Variation of bending strain with pressure at 
the noted locations around the centerline of cylinder 23. 



Cylinder No. 23 
[9O2/OW 




♦ Radial deflections shown at 1 50 psi, magnified 200 times. 

Figure C. 16. . Shape variation due to bending strain for 
Cylinder 23. 
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APPENDIX D 



SPREADSHEETS OF STRAIN GAGE READINGS 
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STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST 



cylinder # layup test date 

6 [903/03]s 2/19/92 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


8 


-2 


3 


-8 


2 


6 


3 


-11 


1 


9 


25 


-80 


86 


-174 


-131 


-145 


-155 


-88 


-144 


-189 


-68 


50 


-124 


29 


-297 


-239 


-257 


-295 


-145 


-236 


-352 


-82 


75 


-136 


0 


-413 


-367 


-343 


-428 


-170 


-314 


-510 


-68 


100 


-160 


-33 


-539 


-445 


-410 


-568 


-181 


-387 


-675 


-33 


no 


-141 


-35 


-575 


-468 


-427 


-634 


-177 


-422 


-784 


2 


120 


-147 


-47 


-631 


-505 


-439 


-690 


-161 


-438 


-829 


37 


130 


-148 


-61 


-695 


-550 


-455 


-767 


-145 


-466 


-925 


81 


140 


-148 


-71 


-753 


-576 


-465 


-852 


-118 


-490 


-1030 


135 


150 


-151 


-91 


-775 


-627 


-472 


-995 


-73 


-527 


-1201 


223 


160 


-149 


-96 


-680 


-645 


-466 


-1056 


-41 


-538 


-1275 


276 


165 


-141 


-95 


-668 


-652 


-452 


-1104 


-14 


-542 


-1336 


321 


170 


-141 


-105 


-701 


-672 


-441 


-1181 


24 


-553 


-1425 


380 


175 


-137 


-no 


-755 


-687 


-427 


-1243 


70 


-557 


-1521 


448 


180 

183 


-132 


-118 


-849 


-711 


-401 


-1338 


115 


-564 


-1613 

-1730 


517 



2. 


Intercept at 


P=0 using 


values 


at 25 and 50 


psi . 










A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


I NT 


-36 


143 


-51 


-23 


-33 


-15 


-31 


-52 


-26 


-54 


3. 


Zeroed 


Circumferential 


Strain 


Data 


V 

(using 


intercept ) . 






Pres A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-44 


-57 


-123 


■108 


■112 


-140 


-57 


-92 


-163 


-14 


50 


-88 


-114 


-246 


216 


224 


-280 


-114 


-184 


-326 


-28 


75 


-100 


-143 


-362 


344 


310 


-413 


-139 


-262 


-484 


-14 


100 


-124 


-176 


-488 


422 


377 


-553 


-150 


-335 


-649 


21 


no 


-105 


-178 


-524 


445 


394 


-619 


-146 


-370 


-758 


56 


120 


-111 


-190 


-580 


482 


406 


-675 


-130 


-386 


-803 


91 


130 


-112 


-204 


-644 


527 


422 


-752 


-114 


-414 


-899 


135 


140 


-112 


-214 


-702 


553 


432 


-837 


-87 


-438 


-1004 


189 


150 


-115 


-234 


-724 


604 


4 39 


-980 


-42 


-475 


-1175 


277 


160 


-113 


-239 


-629 


622 


433 


-1041 


-10 


-486 


-1249 


330 


165 


-105 


-238 


-617 


629 


419 


-1089 


17 


-490 


-1310 


375 


170 


-105 


-248 


-650 


649 


408 


-1166 


55 


-501 


-1399 


434 


175 


-101 


-253 


-704 


664 


394 


-1228 


101 


-505 


-1495 


502 


180 


-96 


-261 


-798 


688 


368 


-1323 


146 


-512 


-1587 


571 



183 -1704 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder H layup test date 

9 (903/03]s 10/31/91 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


-22 


-17 


-18 


3 


1 


-15 


-23 


-18 


-21 


3 


25 


-104 


-98 


-312 


-145 


-184 


-263 


-147 


-156 


-194 


-130 


50 


-137 


-149 


-416 


-259 


-284 


-472 


-213 


-260 


-343 


-232 


75 


-189 


-197 


-444 


-382 


-412 


-673 


-288 


-376 


-496 


-324 


100 


-247 


-256 


-513 


-454 


-526 


-856 


-339 


-471 


-645 


-381 


110 


-245 


-275 


-535 


-470 


-562 


-936 


-344 


-506 


-699 


-390 


120 


-274 


-293 


-574 


-496 


-613 


-1035 


-358 


-547 


-769 


-406 


130 


-286 


-316 


-610 


-517 


-656 


-1134 


-359 


-583 


-836 


-415 


140 


-291 


-333 


-651 


-532 


-689 


-1240 


-357 


-616 


-911 


-419 


150 


-308 


-358 


-714 


-550 


-734 


-1375 


-343 


-659 


-1006 


-421 


160 


-318 


-387 


-788 


-554 


-762 


-1564 


-302 


-708 


-1119 


-402 


170 


-327 


-403 


-868 


-538 


-773 


-1786 


-230 


-748 


-1234 


-350 


180 


-328 


-422 


-1017 


-489 


-757 


-2077 


-104 


-832 


-1392 


-251 


190 

200 


-325 


-443 


-1472 


-265 


-649 


-2743 

-4279 


408 

2152 


-1071 

-2583 


-1653 


65 



2. Intercept at 


P=0 using values 


at 25 


and 50 


psi . 








A* B* 


C 


D 


E 


F 


G 


H 


J 


K 


INT -71 -47 


-208 


-31 


-84 


-54 


-81 


-52 


-45 


-28 



3. Zeroed 


Circumferential 


Strain Data 


(using 


intercept ) . 






Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-33 


-51 


-104 


-114 


-100 


-209 


-66 


-104 


-149 


-102 


50 


-66 


-102 


-208 


-228 


-200 


-418 


-132 


-208 


-298 


-204 


75 


-118 


-150 


-236 


-351 


-328 


-619 


-207 


-324 


-451 


-296 


100 


-176 


-209 


-305 


-423 


-442 


-802 


-258 


-419 


-600 


-353 


no 


-174 


-228 


-327 


-439 


-478 


-882 


-263 


-454 


-654 


-362 


120 


-203 


-246 


-366 


-465 


-529 


-981 


-277 


-495 


-724 


-378 


130 


-215 


-269 


-402 


-486 


-572 


-1080 


-278 


-531 


-791 


-387 


140 


-220 


-286 


-443 


-501 


-605 


-1186 


-276 


-564 


-866 


-391 


150 


-237 


-311 


-506 


-519 


-650 


-1321 


-262 


-607 


-961 


-393 


160 


-247 


-340 


-580 


-523 


-678 


-1510 


-221 


-656 


-1074 


-374 


170 


-256 


-356 


-660 


-507 


-689 


-1732 


-149 


-696 


-1189 


-322 


180 


-257 


-375 


-809 


-458 


-673 


-2023 


-23 


-780 


-1347 


-223 


190 

200 


-254 


-396 


-1264 


-234 


-565 


-2689 

-4225 


489 

2233 


-1019 

-2531 


-1608 


93 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder # layup test date 

12 [904/02]s 11/15/91 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


-5 


2 


17178 


-14 


-6 


1 


-4 


1 


-2 


-4 


25 


-104 


-94 


17194 


-62 


-115 


-94 


-96 


-130 


-71 


-92 


50 


-151 


-142 


17247 


-66 


-211 


-137 


-171 


-395 


-120 


-176 


75 


-203 


-208 


17298 


-69 


-354 


-206 


-269 


-671 


-181 


-284 


100 


-244 


-247 


17346 


-58 


-459 


-148 


-303 


-924 


-202 


-314 


no 


-257 


-250 


17370 


-36 


-497 


-151 


-296 


-1018 


-198 


-314 


120 


-274 


-260 


17390 


-18 


-573 


-170 


-298 


-1188 


-198 


-323 


130 


-289 


-260 


17414 


10 


-661 


-187 


-274 


-1364 


-192 


-324 


140 


-296 


-247 


17434 


37 


-733 


-202 


-242 


-1538 


-174 


-317 


150 


-306 


-225 


17453 


81 


-836 


-214 


-185 


-1783 


-146 


-303 


160 


-323 


-174 


17474 


163 


-1010 


-251 


-42 


-2229 


-71 


-261 


165 


-322 


-123 


17486 


246 


-1104 


-264 


66 


-2475 


7 


-216 


170 


-332 


-32 


17500 


389 


-1334 


-313 


285 


-3012 


152 


-131 


175 


-349 


54 


17511 


543 


-1522 


-377 


530 


-3522 


313 


-45 


180 

181 


-419 


251 


17521 


985 


-1889 


-565 


1182 


-5200 

-5800 


840 


202 



2. 


Intercept at 


P=0 using 


values at 25 


and 50 


psi . 










A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


INT 

3 . 


-57 

Zeroed 


-46 

Circumferential 


-58 -19 

Strain Data 


-51 

(using 


-21 135 

intercept ) . 


-22 


-8 


Pres A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 




0 


0 


0 


0 


0 


0 


0 


25 


-47 


-48 


b 


-4 


-96 


-43 


-75 


-265 


-49 


-84 


50 


-94 


-96 


a 


-8 


-192 


-86 


-150 


-530 


-98 


-168 


75 


-146 


-162 


d 


-11 


-335 


-155 


-248 


-806 


-159 


-276 


100 


-187 


-201 




0 


-440 


-97 


-282 


-1059 


-180 


-306 


no 


-200 


-204 


g 


22 


-478 


-100 


-275 


-1153 


-176 


-306 


120 


-217 


-214 


a 


40 


-554 


-119 


-277 


-1323 


-176 


-315 


130 


-232 


-214 


g 


68 


-642 


-136 


-253 


-1499 


-170 


-316 


140 


-239 


-201 


e 


95 


-714 


-151 


-221 


-1673 


-152 


-309 


150 


-249 


-179 




139 


-817 


-163 


-164 


-1918 


-124 


-295 


160 


-266 


-128 




221 


-991 


-200 


-21 


-2364 


-49 


-253 


165 


-265 


-77 




304 


-1085 


-2 13 


87 


-2610 


29 


-208 


170 


-275 


14 




447 


-1315 


-262 


306 


-3147 


174 


-123 


175 


-292 


100 




601 


-1503 


-326 


551 


-3657 


335 


-37 


180 

181 


-362 


297 


1043 


-1870 


-514 


1203 


-5335 

-5935 


862 


210 



132 . 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder # 
13 



layup 

[904/02]8 



test date 
10/17/91 



1. Raw Circumferential 


Strain 


Gage Data. 










Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


-4 


-6 


1 


-2 


-5 


-1 


-6 


-5 


5 


4 


25 


-120 


-85 


-114 


-102 


-122 


10243 


-85 


-215 


-69 


-109 


50 


-192 


-139 


-186 


-202 


-289 


-182 


-121 


-435 


-116 


-226 


75 


-267 


-174 


-249 


-261 


-366 


1009 


-107 


-633 


-80 


-251 


100 


-306 


-163 


-316 


-309 


-483 


1131 


-80 


-962 


-38 


-304 


no 


-314 


-160 


-357 


-328 


-542 


1100 


-53 


-1123 


-7 


-313 


120 


-325 


-146 


-397 


-340 


-606 


1133 


-11 


-1302 


35 


-318 


130 


-341 


-123 


-451 


-341 


-696 


1133 


48 


-1549 


94 


-332 


140 


-364 


-95 


-516 


-320 


-836 


1167 


146 


-1858 


205 


-348 


150 


-378 


-30 


-595 


-271 


-1024 


1201 


291 


-2298 


367 


-361 


160 


-376 


105 


-732 


-138 


-1311 


1355 


568 


-3010 


659 


-354 


170 

173 


-371 


459 


-1093 


316 


-1947 


1602 


1348 


-4699 

-7500 


1378 


-274 



2. 


Intercept 


at 


P=0 using values at 25 


and 


50 


psi . 










A* 


B* 


C D E 


F 




G 


H 


J 


K 


I NT 


-48 


-31 


-42 -2 45 






-49 


5 


-22 


8 



3. Zeroed 


Circumferential 


Strain Data 


(using 


intercept ) . 






Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 




0 


0 


0 


0 


25 


-72 


-54 


-72 


-100 


-167 


b 


-36 


-220 


-47 


-117 


50 


-144 


-108 


-144 


-200 


-334 


a 


-72 


-440 


-94 


-234 


75 


-219 


-143 


-207 


-259 


-411 


d 


-58 


-638 


-58 


-259 


100 


-258 


-132 


-274 


-307 


-528 




-31 


-967 


-16 


-312 


no 


-266 


-129 


-315 


-326 


-587 


9 


-4 


-1128 


15 


-321 


120 


-277 


-115 


-355 


-338 


-651 


a 


38 


-1307 


57 


-326 


130 


-293 


-92 


-409 


-339 


-741 


9 


97 


-1554 


116 


-340 


140 


-316 


-64 


-474 


-318 


-881 


e 


195 


-1863 


227 


-356 


150 


-330 


1 


-553 


-269 


-1069 




340 


-2303 


389 


-369 


160 


-328 


136 


-690 


-136 


-1356 




617 


-3015 


681 


-362 


170 

173 


-323 


490 


-1051 


318 


-1992 




1397 


-4704 

-7505 


1400 


-282 



/3S 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder # lay\jp test date 

16 [902/(+_30)2]s 11/5/91 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


F* 


B 


C 


D 


E 


G 


H 


J 


K 


0 


1 


4 


10 


-15343 


6 


10 


13 


-8 


-2 


-7 


25 


-23 


-43 


-136 


-15504 


-132 


-139 


-186 


-105 


-184 


-185 


50 


-98 


-33 


-303 


-15574 


-310 


-310 


-385 


-218 


-361 


-347 


75 


-132 


-1 


-428 


-15625 


-445 


-431 


-610 


-325 


-536 


-514 


100 


-161 


24 


-546 


-15676 


-583 


-564 


-861 


-391 


-721 


-687 


no 


-158 


58 


-577 


-15700 


-618 


-606 


-976 


-404 


-794 


-745 


120 


-156 


65 


-589 


-15718 


-646 


-635 


-1097 


-391 


-868 


-815 


130 


-164 


27 


-592 


-15743 


-671 


-662 


-1267 


-376 


-959 


-892 


140 


-166 


-52 


-561 


-15779 


-712 


-675 


-1541 


-322 


-1095 


-1027 


150 


-124 


-124 


-387 


-15825 


-691 


-588 


-1912 


-154 


-1241 


-1187 


160 


7 


-269 


269 


-15915 


-581 


-227 


-2674 


350 


-1449 


-1545 



165 -5200 



2. Intercept at P=0 using values at 25 and 50 psi . 





A* 


F* 


B C ^ 


D 


E 


G 


H 


J 


K 


I NT 


52 


-53 


31 


46 


32 


13 


8 


-7 


-23 



3. Zeroed 


Ci rcumferent ial 


Strain 


1 Data 


(using 


intercept ) . 






Pres 


A* 


F* 


B 


C 


D 


E 


G 


H 


J 


K 


0 


0 


0 


0 




0 


0 


0 


0 


0 


0 


25 


-75 


10 


-167 


b 


-178 


-171 


-199 


-113 


-177 


-162 


50 


-150 


20 


-334 


a 


-356 


-342 


-398 


-226 


-354 


-324 


75 


-184 


52 


-459 


d 


-491 


-463 


-623 


-333 


-529 


-491 


100 


-213 


77 


-577 




-629 


-596 


-874 


-399 


-714 


-664 


no 


-210 


111 


-608 


g 


-664 


-638 


-989 


-412 


-787 


-722 


120 


-208 


118 


-620 


a 


-692 


-667 


-1110 


-399 


-861 


-792 


130 


-216 


80 


-623 


g 


-717 


-694 


-1280 


-384 


-952 


-869 


140 


-218 


1 


-592 


e 


-758 


-707 


-1554 


-330 


-1088 


-1004 


150 


-176 


-71 


-418 




-737 


-620 


-1925 


-162 


-1234 


-1164 


160 

165 

170 

175 

180 

183 


-45 


-216 


238 




-627 


-259 


-2687 

-5213 


342 


-1442 


-1522 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder # layup test date 

17 [902/(4_30)2]s 11/8/91 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


-49 


-146 


-42 


-16 


-15 


8 


-13 


-47 


-31 


-42 


25 


-92 


-192 


-415 


-237 


-132 


-350 


-161 


-557 


-210 


-226 


50 


-131 


-211 


-1322 


-355 


-279 


-476 


-264 


-332 


-327 


-367 


0 


30 


-101 


-642 


-6 


93 


-26 


52 


-22 


31 


-13 


25 


-94 


-170 


-1210 


-226 


-113 


-255 


-143 


-692 


-178 


-203 


50 


-152 


-211 


-1477 


-396 


-257 


-397 


-277 


-323 


-340 


-413 


75 


-162 


-234 


-1643 


-522 


-414 


-582 


-391 


-412 


-484 


-448 


100 


-192 


-262 


-1839 


-650 


-598 


-748 


-518 


-571 


-638 


-543 


no 


-197 


-267 


-1937 


-697 


-698 


-827 


-573 


-641 


-713 


-600 


120 


-208 


-274 


-2009 


-721 


-768 


-880 


-603 


-692 


-759 


-629 


130 


-212 


-277 


-2061 


-752 


-839 


-941 


-643 


-751 


-810 


-665 


140 


-240 


-281 


-2173 


-789 


-974 


-1005 


-706 


-848 


-878 


-726 


150 


-251 


-268 


-2228 


-785 


-1052 


-1026 


-731 


-912 


-902 


-782 


160 


-285 


-256 


-2338 


-763 


-1191 


-1068 


-769 


-1026 


-907 


-839 


170 


-368 


-191 


-2508 


-635 


-1476 


-1077 


-791 


-1280 


-794 


-976 


175 


-418 


-112 


-2589 


-470 


-1708 


-1013 


-774 


-1506 


-602 


-1104 


180 


-686 


227 


-2761 


100 


-2483 


-671 


-704 


-2285 


1062 


-1606 


183 










-2825 












2. Intercept at 


P=0 using values at 25 


and 50 


psi . 










A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


I NT 


-36 


-129 


-943 


-56 


31 


-113 


-9 


-61 


-16 


7 


















(-192025) 




3. Zeroed 


Circumferential 


Strain Data 


( using 


intercept ) . 






Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-56 


-63 


528 


-181 


-163 


-237 


-152 


-496 


-194 


-233 


50 


-95 


-82 


-379 


-299 


-310 


-363 


-255 


-271 


-311 


-374 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-58 


-41 


-267 


-170 


-144 


-142 


-134 


-131 


-162 


-210 


50 


-116 


-82 


-534 


-340 


-288 


-284 


-268 


-262 


-324 


-420 


75 


-126 


-105 


-700 


-466 


-445 


-469 


-382 


-351 


-468 


-455 


100 


-156 


-133 


-896 


-594 


-629 


-635 


-509 


-510 


-622 


-550 


no 


-161 


-138 


-994 


-641 


-729 


-714 


-564 


-580 


-697 


-607 


120 


-172 


-145 


-1066 


-665 


-799 


-767 


-594 


-631 


-743 


-636 


130 


-176 


-148 


-1118 


-696 


-870 


-828 


-634 


-690 


-794 


-672 


140 


-204 


-152 


-1230 


-733 


-1005 . 


-892 


-697 


-787 


-862 


-733 


150 


-215 


-139 


-1285 


-729 


-1083 


-913 


-722 


-851 


-886 


-789 


160 


-249 


-127 


-1395 


-707 


-1222 


-955 


-760 


-965 


-891 


-846 


170 


-332 


-62 


-1565 


-579 


-1507 


-964 


-782 


-1219 


-778 


-983 


175 


-382 


17 


-1646 


-414 


-1739 


-900 


-765 


-1445 


-586 


-nil 


180 


-650 


356 


-1818 


156 


-2514 


-558 


-695 


-2224 


1078 


-1613 



183 -2856 



135 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST 



cy 1 i nder 
18 



layup test date 

[902/(+_30)2]s 12/11/91 



1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


-10 


-6 


4 


2 


3 


2 


-3 


-8 


3 


7 


25 


-61 


-42 


-26 


-106 


-285 


79 


-153 


-155 


-143 


-164 


50 


-82 


-65 


-112 


-160 


-580 


1260 


-306 


-306 


-259 


-296 


75 


-117 


-83 


-259 


-186 


-914 


1100 


-447 


-448 


-348 


-414 


100 


-161 


-94 


-427 


-179 


-1298 


1000 


-588 


-609 


-425 


-536 


no 


-177 


-87 


-509 


-145 


-1485 


993 


-634 


-678 


-433 


-566 


120 


-203 


-78 


-602 


-104 


-1691 


975 


-677 


-762 


-445 


-612 


130 


-241 


-64 


-709 


-32 


-1929 


969 


-711 


-867 


-439 


-661 


140 


-293 


-46 


-841 


69 


-2221 


967 


-733 


-1011 


-416 


-721 


150 


-366 


-15 


-999 


229 


-2576 


976 


-724 


-1199 


-351 


-784 


160 


-485 


71 


-1249 


545 


-3111 


1009 


-631 


-1541 


-191 


-869 


170 


-787 


263 


-1815 


1314 


-4174 


1067 


-269 


-2385 


280 


-976 



175 -5200 



2. Intercept at P=0 using values at 25 and 50 psi . 





A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


INT -40 

3. Zeroed 


-19 60 

Ci rcumf erent ial 


-52 10 

Strain Data 


-1102 
( using 


0 -4 

intercept ) . 


-27 


-32 


Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-21 


-23 


-86 


-54 


-295 


1181 


-153 


-151 


-116 


-132 


50 


-42 


-46 


-172 


-108 


-590 


2362 


-306 


-302 


-232 


-264 


75 


-77 


-64 


-319 


-134 


-924 


2202 


-447 


-444 


-321 


-382 


100 


-121 


-75 


-487 


-127 


-1308 


2102 


-588 


-605 


-398 


-504 


no 


-137 


-68 


-569 


-93 


-1495 


2095 


-634 


-674 


-406 


-534 


120 


-163 


-59 


-662 


-52 


-1701 


2077 


-677 


-758 


-418 


-580 


130 


-201 


-45 


-769 


20 


-1939 


2071 


-711 


-863 


-412 


-629 


140 


-253 


-27 


-901 


121 


-2231 


2069 


-733 


-1007 


-389 


-689 


150 


-326 


4 


-1059 


281 


-2586 


2078 


-724 


-1195 


-324 


-752 


160 


-445 


90 


-1309 


597 


-3121 


2111 


-631 


-1537 


-164 


-837 


170 

175 


-747 


282 


-1875 


1366 


-4184 
-5210 ' 


2169 


-269 


-2381 


307 


-944 



I2>(f 



STRAIN GAGE DATA FROM COMPOSITE CYLINDER HYDROSTATIC TEST. 



cylinder # layup test date 

23 [902/0]4t 11/22/91 

1. Raw Circumferential Strain Gage Data. 



Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


5 


-331 


10 


-4 


28 


26 


11 


-9 


5 


-6 


25 


-96 


-414 


-82 


-92 


-96 


107 


-99 


-89 


-134 


-76 


50 


-147 


-471 


-88 


-218 


-175 


81 


-163 


-119 


-241 


-112 


75 


-191 


-529 


-68 


-321 


-254 


5 


-245 


-144 


-387 


- 142 


100 


-199 


-581 


-104 


-436 


-303 


-77 


-312 


-137 


-601 


-139 


no 


-199 


-601 


-126 


-480 


-299 


-97 


-329 


-106 


-686 


-99 


120 


-211 


-635 


-171 


-565 


-273 


-143 


-346 


-55 


-805 


-33 


130 


-188 


-661 


-203 


-645 


-209 


-233 


-363 


25 


-1037 


76 


140 


-131 


-729 


-332 


-857 


50 


-488 


-350 


263 


-1524 


4 30 


150 

152 


59 


-926 


-834 


-1306 


1028 


-1305 


-107 


849 


-2559 

-4400 


1362 



2. Intercept at 


P=0 using values 


at 25 


and 50 


psl . 








A* B* 


C 


D 


E 


F 


G 


H 


J 


K 


INT -45 -357 


-76 


34 


-17 


133 


-35 


-59 


-27 


-40 



3. Zeroed 


Cl rcumf erent 1 al 


Strain 


1 Data 


(using 


intercept ) . 






Pres 


A* 


B* 


C 


D 


E 


F 


G 


H 


J 


K 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


25 


-51 


-57 


-6 


-126 


-79 


-26 


-64 


-30 


-107 


-36 


50 


-102 


-114 


-12 


-252 


-158 


-52 


-128 


-60 


-214 


-72 


75 


-146 


-172 


8 


-355 


-237 


-128 


-210 


-85 


-360 


-102 


100 


-154 


-224 


-28 


-470 


-286 


-210 


-277 


-78 


-574 


-99 


no 


-154 


-244 


-50 


-514 


-282 


-230 


-294 


-47 


-659 


-59 


120 


-166 


-278 


-95 


-599 


-256 


-276 


-311 


4 


-778 


7 


130 


-143 


-304 


-127 


-679 


-192 


-366 


-328 


84 


-1010 


116 


140 


-86 


-372 


-256 


-891 


67 


-621 


-315 


322 


-1497 


470 


150 

152 


104 


-569 


-758 - 


1340 


1045 


-1438 


-72 


908 


-2532 

-4373 


1402 
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APPENDIX E 



ABD STIFFNESS MATRICES 



138 



K 

((90)4/(0)2]S| r - 2.755", 1-18 



Buckling Pressure for a cylinder with dlicenslons ofi 
r - 2.75, 1 - 18.00, t « 0,0720 -> I - 1633. , r/t - 38. 

and a load ratio of: Ky - 2.00 Nx - 1.00 



Per - 


218.2 


NX - 


300.6 M - 


1 N 


- 2 


ABAQUS Per - 


242.5 












A 


BO Matrix: 
















5.3129E+05 


3.0428E4^04 


0.0000 




0.0000 


0.0000 


0.0000 




3.0428E+04 


9.5390E+05 


0.0000 




0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


7.2000E+04 


0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


0.0000 




1 67.23 


13.15 


0.0000 




0.0000 


0.0000 


0,0000 




13.15 


574.4 


0.0000 




0.0000 


0.0000 


0.0000 




1 0.0000 


0.0000 


31.10 


I(90)3/(0)3)s; r 


- 2.755", 1 


- 18" 










Buckling Pressure for a 


cylinder 


with dlaenslons 


of: 




r - 2.75, 


1 - 18.00, 


t - 0.0720 


-> Z - 1633. , 


r/t - 38. 


and a 


load ratio of: Ky - 


2.00 NX - 


1.00 




Per - 


236,6 


NX - 


325.9 M - 


1 N 


- 2 


ABAQUS Per - 


230.5 












ABD Matrix: 
















7.4260E+05 


3.0428E+04 


0.0000 




0.0000 


0.0000 


0.0000 




3.0428E^04 


7.4260E+05 


0.0000 




0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


7.2000E+04 


0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


0.0000 




115.4 


13.15 


0.0000 




0.0000 


0.0000 


0.0000 




13.15 


526.2 


0.0000 




0.0000 


0.0000 


0.0000 




0.0000 


0.0000 


31.10 



((90)2/(±30)2]s; r - 2.755", 1 - 18" 

Buckling Pressure for a cylinder with dl»enslons of: 
r - 2.75, 1 - 18.00, t - 0.0720 -> Z - 1633. , r/t - 38. 

and a load ratio of: Ny - 2.00 “ Nx - 1.00 

Per - 206.6 Nx - 284.6 M - 1 N - 3 

ABAQUS Per - 205.4 

ABD Matrix: 



6.0055E+05 

1.7247E+05 

0.0000 


1.7247E+05 

6.0055E+05 

0.0000 


0.0000 

0.0000 

2.1404E+05 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


0.0000 

0.0000 

0.0000 


141.4 

40.42 

19.08 


40.42 

445.7 

7.271 


19.08 

7.271 

58.38 



r/t - 38 



Buckling Pressura for a cylinder with dlncnslons oft 
r - 2.75, 1 - 18.00, t - 0.0720 -> Z - 1633. 

and a load ratio oft Ny - 2.00 Kx - 1.00 

Per - 184.1 Hx - 253.7 M - 1 N - 3 

ABAQUS Per - 178,8 



ABD Matrix: 



5 . 3129 E +05 


3 . 0428 E +04 


0.0000 


1 2536 . 


0.0000 


0.0000 


3 . 0428 E '(-04 


9 . 5390 E +05 


0.0000 


0.0000 


- 2536 . 


0.0000 


0.0000 


0.0000 


7 . 2000 E +04 1 


1 0,0000 


0.0000 


0.0000 


2536 . 


0.0000 


0.0000 


234.6 


13.15 


0.0000 


0.0000 


- 2536 . 


0.0000 


13.15 


407.0 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


31.10 



Jones (90,90,90,60,30,0,0,15,30,60,75,90); r - 2.755", 1 - 18" 
Buckling Pressure for a cylinder with dimensions of: 



r - 2.75, 


1 - 18.00 


, t - 0.0720 


-> z - 


1633 . 


, r/t - 38. 


and a 


load ratio 


Of: Hy - 


2.00 Kx - 


1.00 




Per - 


235.9 


NX - 324.9 M 


-IN 


- 3 


ABAQUS Per « 


213.4 










ABD Matrix: 














S. 540 eE-f 05 


1 . 1329 E +05 


1 . 1791 E +05 


1286 . 


426.1 


565.4 




1 . 1329 E +05 


7 . 6539 E +05 


1 . 1791 E +05 


426.1 


- 2138 . 


934.5 




1 . 1791 E +05 


1 . 1791 E 405 


1 . 5486 E +05 


565.4 


934.5 


426.1 




1286 . 


426.1 


565.4 


83.12 


35.45 


22.95 




426.1 


- 2138 . 


934.5 


35.45 


513.9 


43.62 




565.4 


934.5 


426.1 


22.95 


43.62 


53.40 




Thesis 
M5933 
c. 1 



Miller 



Hydrostatic testing 
analysis of graphite! 
cylinders . 



and 

epoxy 



Th(^sis 

M5933 Miller 

c.l Hydrostatic testing and 

analysis of graphite-epo5^ 
cylinders. 



